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^ . Abstract 

>: 

Using the expansions of the heavy meson decay widths in the heavy quark mass 

^T) ■ and QCD sum rules for estimates of corresponding matrix elements, we calculate the 

jj±,o,s decay widths and the B ±,0,s lifetime differences. The results for D mesons are 

in agreement with the data, while it is predicted that [T(B°) — T(B~)]/T B ~ 4%, 

and the lifetime difference of the B° and B s mesons is even smaller. The role of the 

& weak annihilation and Pauli interference contributions to the lifetime differences are 

1 ■ described in detail. In the course of self-consistent calculations the values of many 

*^h. parameters crucial for calculations with charmed and beauty mesons are found. In 

particular, the perturbative pole quark masses are: M c ~ 1.65 GeV, Mb ~ 5.04 GeV , 

j> \ and the decay constants are: f D (M c ) ~ 165 MeV , /b(M 6 ) ~ 113 MeV . It is also 

shown that the nonfactorizable corrections to the B° — B° mixing are large, Bb — (1 — 

18%) . The values of the unitarity triangle parameters are found which are consistent 
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X 



a 



with these results and the data available (except for the NA31 result for the e'/e which 
is too large): A ~ 0.22 , A ~ 0.825 , p ~ -0.4 , 77 ~ 0.2 . 
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1 . Introduction 

It is a long standing challenge for theory to calculate the D and B meson decay 
widths. On the qualitative side, two mechanisms were invoked to explain the pattern 
of the D meson lifetime differences: weak annihilation (WA) [l|], [FJ, ||, and Pauli in- 
terference (PI) [U , || , , • As for WA, it was expected that because an admixture 
of the wave function component with an additional gluon or the emission of a pertur- 
bative gluon, both remove a suppression due to helicity conservation (which leads to 
Br (it — > ev)jBr (n — > \xv) ~ 10 -4 ), the D° meson decay width is enhanced. On the 
other hand, it was expected that the destructive Pauli interference of two d-quarks 
(the spectator and those from a final state) suppresses the D^ meson decay width. As 
for WA, there were no reliable calculations at all. For PI, simple minded estimates 
(see sect. 8) give too large an effect which results in a negative D^ decay width. 

For B mesons, it was clear qualitatively that all the above effects which are of 
a pre-asymptotic nature and die off at Mq — > oo, will be less important. However, 
because the patern of the D meson lifetime differences was not really explained and 
well understood, this prevented to obtain reliable estimates of the B meson lifetime 
differences, and only order of magnitude estimates are really available: [5T(B)/Tb] '■ 
[5T(D)/T D ] ~ OUlMt/flMl) ~ O(H)- 1 ). 

Moreover, as for WA contributions through perturbative gluon emission (which is 
formally a leading correction ~ 0(Aq C d/Mq) to the deacay width and was expected 
before to be potentially the most important), it has been emphasized recently |§ 
that such contributions are of no help at all because, being large (at least formally 
at Mq — > oo) term by term, they cancel completely in the inclusive widths, both 
0(1/Mq) and 0(1/Mq) terms. [] It will be shown below (see sect. 9) that, nevertheless, 
there are important WA contributions but on the nonperturbative level. 

Considerable progress has been achieved recently in applications of the operator 
product expansion to the calculation of the heavy meson decay width. In particular, 
it was shown that there are no 0(1/ Mq) corrections to the Born term and first non- 
perturbative corrections 0(1/Mq) were calculated explicitly P[, [lQfl. However, these 
contributions are all nonvalence (i.e. one and the same for all D ±,0,s mesons), and so 
have nothing to do with lifetime differences. They are important however for the cal- 



culation of the absolute decay rates. Recently a number of papers appeared [O], [12 



13| , [|n|1, where these results were applied to determine the values of the quark masses, 



M c and M b , and [V^,). 

The purpose of the present work was to calculate the D and B meson decay widths, 
with the main emphasis on the calculation of lifetime differences. It is shown that the 
£)±,o,s lifetimes can be calculated with a reasonable accuracy, and concrete predictions 



1 The separate perturbative WA contributions into the total width have the form: ST /T 'Born ~ 
/|)/ e o > fh ~ Po/Mq, where e Q is the binding energy of the spectator quark in the meson. So, they 
are highly infrared sensitive and singular at e — > 0. It seems clear that there can not be power infrared 
singularities in the inclusive decay width, so that it is not so surprising that all such terms cansel (i.e. 
the terms ~ 0(1/ 1"^) and ~ 0(l/e ), but not loge which describes hybrid log. ) . 



for the B ±,0yS meson lifetime differences are given also. 

The scope of the paper is as follows: 1) Introduction 
2) Definition of the heavy quark mass 
3a) General formulae: c-quark 
3b) General formulae: b-quark 

4) D —>■ e v + X. Determination of M c 

5) Mass formulae. Determination of Mb and A 

6) B — > e v + X. Determination of \V c b\ 

7) Calculation of fo, Jb 

8) Difficulties with naive estimates 

9) Nonfactorizable contributions: gluon condensates 

10) Nonfactorizable contributions: quark condensates 

11) Corrections to semileptonic widths 

12) A and A 

13a) Calculation of D ±,0)S decay widths 
13b) Calculation of B ±,0,s lifetime differences 

14) B° - B° mixing 

15) The unitarity triangle 

16) Summary and conclusions 

2. Definition of the heavy quark mass 

In what follows, expansions in powers of the heavy quark mass and some formulae 
of the HQET (Heavy Quark Effective Theory) (see, i.e. the reviews ||i"5| , fi~6H , [|l7l ) 



are used. In the standard approach, the perturbative pole mass, M p , is chosen as the 
HQET expansion parameter. The perturbative pole mass is, clearly, a distinguished 
parameter because it is scheme and gauge independent. So, it is natural that all 
calculations of observable quantities are expressed usually through M p . 

It was pointed out recently |I8||, fll9"|, p0| that the perturbative pole mass is an 
ill-defined quantity and contains an intrinsic uncertainty ~ O(Aqcd), because the 
perturbation theory series for it diverges due to renormalon effects. 

Analogous renormalon divergences are well known and originate from the fact that 
there is an internal scale in the theory: Aqcd~ the position of the coupling constant 
infrared pole. So, this scale reveals itself when perturbative loops are integrated up to 
ki — > (as in the case with the pole mass). But because Aq C d can not appear explicitly 
at any order of renormalized perturbation theory, the perturbative series diverges. 
These infrared singularities are cured usually by nonperturbative contributions, and it 
seems most people are not worrying especially about them. 

To avoid from the beginning the infrared region contributions, it was proposed in 
|i~8|l to use some infrared safe mass parameter, M M , instead of M p . As for a concrete 
choice of M M , two variants are considered usually. The first one (in a spirit of the 
Wilson operator expansion) is to calculate all loops by cutting out the ki < /i region 
contributions. A deficiency of this variant is that it is practically impossible to per- 



form such calculations. The second one is to use the current quark mass normalized 
far off mass shell which suppresses the infrared region contributions. Because this 
current mass differs from M p already at the first loop level, this variant spoils finally 
all the usefullness of the standard HQET approach, replacing the standard decompo- 
sition in powers of 1/M p by the usual perturbation theory decomposition in powers 
of 1/logM^. Besides, because the current quark mass is renormalization scheme and 
gauge dependent, all this is highly inconvenient. 

In connection with this, we describe below another redefinition of M p and introduce 
the "hard pole mass", M , which is free of renormalon singularities by definition and 
differs from M p by only the term ~ ^qcd (not by ~ a s M p ). We emphasize that, 
as far as we do not calculate explicitly high order perturbation theory corrections to 
observable quantities, the whole problem of the renormalon is of abstract interest only 
and can be formally solved on a "hand waving level", i.e. by some redefinitions only, 
without making any real changes in all the formulae available. 

The above "hard pole mass", M , can be connected with a definition of the matrix 
element of the local operator. Although it is clear beforehand, let us demonstrate 



explicitly, using the example considered in [18], that the relevant operator to use for a 



redefinition of M p is Hn g ht, i.e. the light degrees of freedom part of the Hamiltonian. 
With this purpose, let us consider the Hamiltonian and take its average over the heavy 
meson state, \Mh >'■ 

< M H \ H tot \M H >= M^<M H \ (Q Q)„ \M H > + 

< M H \ ^(E 2 + H 2 )^\M H > + <H q > +0(^ 2 /M) , (1) 

where \i serves as an upper cut off in the matrix elements of operators and as lower cut 
off in Mfj,, H q is the light quark Hamiltonian, < M H \(QQ) fl \M H >= 1 + 0(fi 2 /M 2 ) . 
In the example considered in |TH] (in the static limit): 
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«*>=i^ *■+*„■ -»■-* ■ (2) 

and when a s in Eq.(2) is substituted by running a s (k 2 ) in the form: 

a,(?) = a,tf)± {^ll 0g Q\ (3) 

this leads to a divergent perturbative series: 

OO L 

5M = nY,as n C n , C n ~{^) n n\. (4) 



On the other hand, it is not difficult to see that the insertion of the vertex E 2 /2 into 
the same diagram which gives SM, adds the contribution: 

, , .1- 9| „, . 167T [ d A k , ... / — k 2 \ ritr .. 

< M H \ -& \M H >= -,_ / — a, m [^- j = SM , (5) 

which cancels the renormalon contribution in Eq.(l ) and leaves the pole mass M p 
instead of M M . 

Let us proceed now with a standard on mass shell renormalization scheme and 
rewrite Eq.(l ) in the form: 

< M H \ H tot \M H >= M p < M H \QQ\M H > + < M H \ H Hght \M H > + 

^ + 0(A^/M 2 ), (6) 

where @ (see sect. 5 for the explicit form of 5fh 2 ): 

<M H \QQ\M H >=l-5rh 2 /M 2 + 0(A^ CD /M 3 ), 6rh 2 = 0(A 2 QCD ) . (7) 

It is not difficult to check that: 

< M H \ H light \M H >= A (l + 0(A 2 QCD /M 2 )) , (8) 

where A is a finite number ~ Aqcd, independent of M. Clearly, it can be understood 
as the light quark self energy plus the binding energy. 

Now, let us redefine the quark pole mass, M p , and introduce the "hard pole mass" , 
M : 

M p = M (1 + Ax + A 2 + ■ • •) , (9) 

where A x is defined formally, for instance, as the leading asymptotic part of the 
divergent perturbative series for M p : 

^=f<(M )CM=A^, C^-^Jnl, (10) 

and A\ ~ (1 + 0(a s )), Aq ~ l/a s (M ). A 2 in Eq.(9) is defined analogously through 
a subleading part of the series, so that A 2 = A 2 Aq CD /M 2 , etc. So, by definition, M 
in Eq.(9) is free of divergent parts of renormalon contributions. What is important, is 
that the "bad" quantity M p and the "good" quantity M differ not by ~ a s (M )M 
terms, but terms ~ Aq C d only 

2 This is natural because the whole answer has to be independent of /z, and considering formally 
the case when there are no nonperturbative contributions and the matrix element is taken over the 
on mass shell heavy quark state, the total answer is M p . 

3 More precisely, if any meaning can be given to the divergent sum in Eq.(10) , it can also give 
higher order terms ~ Aq CD /M^ + ..., but these can be included into a redifmition of A2, etc. 



Clearly, we can absorb the term Ai into a redefinition of A , the term A2 into a 
redefinition of 5m 2 , etc., so that Eq.(6) can be rewritten as: 

2 5m 2 
< M H \ H tot \M H >= M < M H \QQ\M H > +A + -j f - + 0(A 3 QCD /M 2 ) , (11) 



A = A + A!M , 5m 2 = 5m 2 + A 2 M 2 . (12) 

In connection with the above, let us point out the following. Strictly speaking, 
because the real expansion parameter is Aqcd/M, taking account of higher order 
terms of the perturbative series in a s (M) is not justified without simultaneously tak- 
ing account of power corrections in Aqcd/M. So, even if we have a convergent series: 
J2 a™{M)B ni the term af 1 (M) B Nl becomes 0(Aq C d/M) starting with N± (depending 
on the behaviour of B n at large n); starting with N 2 the term a^ 2 (M) B^ 2 becomes 
0(A 2 QCD /M 2 ), etc. For instance, for B n ~ 1 at large n: JVi ~ [a s (M) logl/a^Af)] -1 , 
while for B n ~ n! : N\ ~ [a g (M)] -1 , which is only slightly larger than those for 
the convergent series. Therefore, in any case, the account of distant terms of the 
perturbative series has to be made only simultaneously with account of power correc- 
tions, and we can always transfer these distant perturbative terms into a redefinition 
of power corrections (remembering especially that we are unable to calculate directly 
these nonperturbative corrections at present). 

Let us turn now to the heavy meson decay width which can be represented in the 
form PI: 



M 5 
p 



F(a s (M p ))+0(A 2 QCD /M 2 )} , (13) 



where F(a s ) represents a perturbative series, and power corrections in Eq.(13) are 
explicitly expressed through the matrix elements of higher dimension operators (see 
sect. 3). 

As was emphasized in |I5[ , a new element in the example considered is that the 
renormalon effects lead to a contribution ~ Aqcd in the pole mass, while there are 
no nonperturbative corrections to V at this level. Correspondingly, the perturbative 
series for F(a s ) is also divergent, and only a product of them in Eq.(13) contains no 
~ Aqcd/M renormalon contributions. 

Now, let us reexpress the pole mass in Eq.(13) through the formally defined "hard 
pole mass", M . Then in the product: (1 + 5 Ai + 0(A 2 )) F(a s (M )) (where Ai is 
understood as the divergent series Eq.(lO)), the leading divergences cancel, so that 
there is no Aqcd/M correction. As for subleading divergencies which give corrections 
~ Aq CD /M 2 , they will be cured in the usual way by explicit nonperturbative matrix 
elements in Eq.(.13) , etc. 

The final result of all the above manipulations can be formulated as follows. 
1. We can use the standard perturbative pole mass, M p , in the HQET to obtain usual 
expansions in powers of 1/M p , putting no attention at the time that it is not a "good" 



quantity. 

2. When dealing with formulae for directly observable quantities, we can reexpress 
M p through the "hard pole mass", M , which is a "good" quantity. All renormalon 
divergencies will be either explicitly canceled in perturbative expansions, or absorbed 
by redefinitions of the original nonperturbative matrix elements. The net effect will 
be that all the original formulae of the 1/M expansions will formally stay intact after 
all the above redefinitions. Moreover, if we calculate explicitly only a few lowest order 
loop corrections, which is the case usually, we need not even explicit expressions for 
the divergent tails of perturbative expansions. It is sufficient to have in mind that all 
the divergent contributions are cured by redefinitions. 

3. What is really important, is that (having in mind all the above redefinitions) we 
can use now our formulae with only a few first loop corrections explicitly calculated for 
a comparison with the data, not worrying much that there may be numerically large 
corrections from distant terms of the perturbative series. 

3a. General formulae: c-quark 

The weak Hamiltonian used in what follows has the form ( T^ = 7^(1 + 75) ): 



H c w (fi = M c 
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V2 



V cs V: d {C 1 (M c )sT^c-uT IM d + C 2 (M c )uT fM c-sT IM d} + --iU) 



where the dots denote corresponding Cabibbo-suppressed contributions. The coeffi- 
cients C± = (Ci ± C 2 )/2 are determined by: 



C ± {^)=C^)ll 



a s (/i 2 ) -a s (M^) 



TX 



P± 



(15) 



C 1 



' a s (/i 2 ) 

M M w) 



4/6 



C^ 2 ) = 1/y/CW) 



(16) 



and the values of p± can be found in |21| . We use below in our calculations: 
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(153-19n/), £ = log(Q 2 //j 2 ) . {It 



One obtains from Eq.(17) (Z = l/a a ): 

b 



Z(^) = Z(Q< 



An 



\og{Q l /n l ) - A log 



'Z(Q 2 



A^ 



Z(/i 2 ) + A 



(19) 



We use as an input the value: 



oc s {M\ 



2 ^ 

W i 



0.118, 



(20) 



which corresponds to: A-^- ~ 200 MeV, A4J= ~ 300 MeV . One obtains then from 



MS 



MS 



Eqs.(19),(20) 0: 

a s (fi = M b = 5.04 GeV) ~ 0.204, a s (/i = M c = 1.65 GeV) ~ 0.310 , (21) 
and the coefficients C{ are equal to: 

C_(M C )~ 1.770, C + (M C ) ~ 0.762 , C X (M C ) ~ 1.266, C 2 (M C ) ~ -0.504.(22) 
The decay width of a hadron containing a heavy quark: 

T= 2W I dX < M H\H w (x)H w (0)\M H >=^ r <H\L eff (0)\H>, (23) 

can be represented in the form of the operator expansion || : 

L e// = £CiOi(0), (24) 

where the leading term is the dimension 3 operator QQ which describes (at Mq — > oo) 
the free quark decay, fig.l. The next term is the dimension 5 operator QaGQ which 
describes an emission of a soft gluon, fig. 2, etc. 



{Fig.l} 

Fig. 1 The Born contribution 

{Fig.2} 

Fig. 2 The soft gluon emission giving rise the correction 0(1 /Mq) 

If we confine ourselves temporarily to these terms only, then all D°' ±,s -mesons will 
have equal decay widths, and (on account of the radiative correction [[H|]) it can be 
represented in the form ||: 



nl — j. ^ 



-L nZn 



'2Ci + C 2 _\ <D\cc\D> 



2M 



D 



,„ j(1 _ 2 | Ag ) +4 |(J|^|W (25) 



Here: 



/, 



2Cl 



rad 



2C\ + Cl 

c 2 _ \ 



I -^7ra s (M c 2 ) + 



2 
3 



12 71 



+ 



2CI + C1 



1 - -Tra a {M c ) + 



(26) 



^The quark masses are found below. 



2 

z = (l-8x + 8x 3 -x 4 - 12x 2 logs) ,Z! = (1- x) 4 , z 2 = (1 - x) 3 , x = -r£,(27) 

GlM 5 < Dice ID > / 1 A \ 

r ^li^' 2M D ^(^- 2 *c-A K ), (28) 

1 <D|c(M c 2 -^)c|D> ^ 1 <D\c^g s a, u G^^c\D> 
G M? <D\cc\D> M? <D\cc\D> ' l ' 



J_ <D\c(jf)c\D> = < (p*) c > 
K 2M C 2 <£>|cc|£>> " 2M C 2 ' l J 

where p c is the c-quark 4-momentum operator, p is its 3- momentum operator (in 
the D-meson rest frame). Let us recall that the radiative corrections in Eq.(26) are 
calculated at m s = 0. 

The lifetime differences appear first on account of the 4-quark operators, and we 

(c) 

°ff 



present here the contributions to SL^L from figs. 3 a, b, c-diagrams (see also Appendix): 



5L tff( M C ) = q ^^{9^ 2 LU + T^( L ; u + L ; u )} +-.., (31) 



L d llu = {S d {cT,d){dT u c) + d {cY,-d){dT u -c) \ , (32) 



\a \a 

d) ( dT v 

2 M 2 



Ll v =ls:(cT,u)(uT v c) + 0:{cT^u){uT^c)\, (33) 



L% = {s o s {cT,s){-sY u c) + 0:{cT^s){-sT^c)\. (34) 



T, 



„„ 3 



-(A M A„- X 2 g^ u ) , (35) 



where A (or A ) is the total 4-momentum of the integrated quark pair. It can be read 
off from each diagram and differ from p c by the spectator quark momenta. 



{Figs. 3a, 3b, 3c} 
Fig. 3 The diagrams contributing into the 4-fermion operators 

8 



The coefficients Sf and Of are: 



no 



3 (c? + cl) 



2C 1 C 2 , 0° d = 2{Cl + C 2 2 ) 



(36) 



SS = 3(C 2 + ic 1 ) 2 , QS = 2C?, 5 s ° = 3(C 1 + ^ 2 ) 2 , O s ° = 2Cf, (37) 
while Ci, C*2 are given in Eq.(22). Therefore: 



Co 


o° d \ 
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3.71 
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0.02 


3.20 
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0°s J ,_ 


-u \ 


v 3.61 


0.51 



The leading contribution at M c — > oo(A ~ A ~ p c ~ M c ) in Eqs.(31)-(34) coincides with 
those obtained before in f|, ||, @. As will be shown below, however, it is important 
for calculations with charmed mesons to account for nonzero momenta of the spectator 
quarks : A ^ A ^ M c . 

The operators entering SL^L are normalized at the point jj, = M c . So, it is reason- 
able to separate out explicitly the "hybrid" log effects |7|1,[p2|,||23l] by renormalizing to 
the point jjl = fi , before the calculation of matrix elements. The 4-quark operators 
in 5L e ff are renormalized as [fj2] (we neglect possible changes in the renormalization 
formulae originating from the presense of the spectator quark momentum operators in 
A): 



Wfiv, 



H=Mq 



t co — 1 2 

- ) bfj, u — - [T co — 1) O^y 



(3? 



Ai=Mo 



(c^V 



Mq 



:i + 



7V, 
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1)V 
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(39) 



s^ = (cT, q )(qr u c), oi v = (sr£q)(qr£c] 



(40) 
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M M Z)J 



(41) 



As for a concrete value of /x , let us point out that the mean value of the vacuum 
quark 4-momentum squared is 



, < 01 qD 2 ,q\0 > < 01 q\ q s a uv Gf lu ^- q 10 > „ , , 

< -k 2 > = A, -J, = — nl - £ — - °- 4 GeV ■ ( 42 ) 



< qq\0 > 



< gg > 



This number determines the characteristic scale of nonperturbative interactions, so 
that the scale /i must not be chosen below of this value. In what follows we will use 
the value: 



^ = o.5 Gel/ 2 , a s (fi 2 = 0.5 GeV 2 ) = 0.5£ 



(43) 



a s (M? 



('■•) 



1.87, r co = rjll 2 ~ 1.37 



Now, after renormalization $L e h has the form: 



where: 



5L { ^ f (j2 = AL(fj, ) + L^yifio) , 



dp 



A!' 



AL{^) = -^{T ltv L% + A i + \'} 



\V U d\ g^v^ L fiU + \V cd \ T^ LV yS s S llu J rO s 0^ u j 



(44) 



(45) 
(46) 

(47) 



A s — \V CS \ T IXV L S IXV + \V US \ g^ v \ yS d S s ^ v + O d O s ^ v 
and the coefficients are: 



Sd o d 
s„ o. 



M=Mo 



-0.73 2.93 
t,.„\ -0.11 2.43 
3.48 -0.26 



-1.07 4.02 
-0.15 3.33 
4.76 -0.36 



(48 



The "penguin nonvalence" (PNV) term, L PNV , originates from the contribution of 
the diagram in fig. 4. 

{Fig. 4} 

Fig. 4 The nonvalence factorizable penguin contribution (PNV) 

With logarithmic accuracy Lp^y can be easily obtained from Eqs.(31)-(34) by 
extracting the term ~ D^G^ U = —g s J u from the light quark operator [qq\. For instance: 



dd] 



1 . A a 

'2 (7 'T 



A a 
dlp—d 



j(7py)^(3) 



1 - nz? /bo , A c 



'" '7„ T ) J?(3)049) 



■^(3) = («7p-2-« + (d) + ( S ; 



(50) 
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where the last step in Eq.(49) describes the evolution of J a from M c to \x . Proceeding 
in this way (and using A ~ A ~ p c , A 2 ~ A 2 ~ M 2 , which is justified for these 
contributions, see fig. 4) we obtain: 



^v(^ )^(i-?u 2/9 )m 2 j;(3) 



J PNV\ 



2vr 



A a 

c7 P y (A^ + iY a75 )c 



(51) 



Ho 



N v = -(A u + A s -2A d )~2.54, 



N n = -(- A, + -A.-2A 



3 V3 



d) -1-62, 



(52) 
(53) 



A i = (S i --O i ). 

6 



(54) 



For the semileptonic decays, SL^ can be easily obtained from Eqs.(45)-(48), (51)) 
by corresponding replacements and has the form : 



crlept _ a j lept i t lept 
0J ^eff = L ^ Ll "T ^PNV i 



(55) 



AL lept 



Gj 

2tt 



T; 



/i V 



X" 



x a 



\V CS \ 2 1.327 cr„5-sr v c- 0.245 cT.— s-sT u —c + 



H c c J- j, ^ u.^t:^ ^j.^ ~~"9 



|Kd| 1.327 cr M d • dr v c- 0.245 cr^—d-dr^—c 
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(56) 
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^ a (3) • c lp (l + i 75 ) 



A" 



II o 



(57) 



3b. General formulae: b-quark 



The relevant part of the weak Hamiltonian has the form: 

H b w (fi = M b ) = ^= V cb V: d { d(M b ) ciy, ■ d~i> + C 2 (M b ) dT ,b ■ cT .u) 



{he 



a 9 (n = M b = 5.04 GeV) ~ 0.204, 



(59) 



C_(M 6 ) ~ 1.385 , C+(M 6 ) ~ 0.855 Ci(M 6 ) ~ 1.120 , C 2 {M b ) ~ -0.265 . (60) 
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For the beauty-hadrons one can obtain SL b e r, from SL^rr by evident replacements: 



SL b eff (M b ) 



2 IT/ 12 



G F \V C b 
2~7" 



{(L u + L d + L c ) + (L C + L 8 + L C )} =m +■■■, (61) 



M, 



X" 



A" 



L u = g^\ 2 1-^H {S° d (bT,u)(uT l/ b) + O d (bT,-u)(uT l/ -b) ,(62) 



A 



L c = g fMU X 1 



2/1 ^ ^(&r,c)(cr^) + o°(6r^c)(cr^fe) ,(63) 



^ = ^(A)|^(&r M rf)(Jr^) + o:(6r M yrf)(Jr,y6)| 



(64) 



^=^(A)|^(6r, s )(sr^)+o:(6r^ s )(sr^6)} 



(65) 



Lc = T^(\)ls s CbT^c)(uT I/ c)+0 s CbT^c)(cT u '—b] 



T 



(66) 



Z c =T#(A)js?(5r A 



c)(cr I/ 6) + oj(5r /i yc)(cr I/ y6; 



(67) 



where the dots in the above formulae denote the Cabibbo-suppressed contributions. 
Here (x = M c 2 /A 2 ) : 



Tni,{\) — - [Xfj, X u — A g ill 



(61 



^(A) = 7,(1 - *) 2 [(>„ A, - X 2 g, u ) (1 + 2x) + - ffp A 2 



;r 



(69) 



T<%\\) = l -VT^Ax~ [(1 + 2x) A„ A„ - (1 - x) X 2 g^ 



The coefficients in Eqs.(61-(67)): 



(70) 



Co 
°d 


o° d \ 


1 


' -0.15 


2.65 


Co 


0: 


* 


0.03 


2.51 


s° 


0° 


., I 


3.20 


0.14 



s s 



H=M, 
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become after the renormalization, Mb — > fi : 

i) for the terms Ld, L s , L u , for which the renormalization factor is 



27/50 i -71 r 

Tbo = % c T^ ~ 1.715, 7? 6c 



« S (M C 2 ) 
a s (M 6 2 ) 



1.52. 



(71) 



Sd 


o d 


&u 


o u 


s s 


o s 



M=Mo 



-0.31 1.71 

r /)( ,j -0.12 1.57 
3.04 -0.80 



-0.53 2.93 
-0.21 2.69 
5.21 -1.37 



ii) for the terms L c , L c , L c , for which the renormalization stops at \x — M c and the 
renormalization factor is: 



n c = Vb 7 J 5 ° - 1-253 : 



(72) 



Sd Od 
5, O, 



ii=M c 



The 



(b) 



-0.23 2.20 
3.12 -0.32 



-0.29 2.75 
3.91 -0.40 



The corresponding term Lp NV has the form: 



L ^ n/ ^3^l M ^r(3) 



J PNV 



2tt 



b^- lp (A + B l5 ) b 



+ J a p (c) 



b^- lp (C + D l5 ) b 



(?3) 



A = X, (A° - A° d ) + A 2 (A° u - A° d ) ~ 0.40; 



(74) 



B = Xi (\a° s - A d ^j + A 2 (±A° U - A d 



0.27, (75) 



C = X 3 (A° + A° u - 2 A° d ) ~ 0.31 ; D = X 3 l-A° + -A° u - 2 A d ~ 0.17 , (76) 



, _ -8/25 9 / i I q -8/25 

A 3 = LJ |^t 6c , A, = A 3 ^ 9 r co , A 2 = ± 1 - 1 + 7 C C 2/9 ] r b0 , (77) 



V =[S°- 1 -0°), J?(3) 



_A a 



(d) 



■W 



6 " V ' p 
For the leptonic decays, B — ► / + z/ + X, / = e, /x, we have: 



_A a 
CyTpC 



«#/ 



AL Zept + L 



£ , r lept 



PNV ? 



(78 



(79) 



AL lept 



n 2 it/ i 2 

Up \V rh\ 



F \ v cb \ 

2vr 



-T,„(X) 



A" 



A a , 



1.225 6T„c • cY v b - 0.169 6T,,— c • cT v — b 



fJn 



?0) 
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T le V l ~ 
^PNV — 



Gl\V cb \ 2 1-% 



-8/25 



2vr 



TboVt 



5 9 Ml 



X 



67,[l + 57 5 )y6 



J?(3) + r- 1 ^ /9 J?( 



/^2 It/ 12 



2vr 



0.046 Ml 



1 \ A a 



J'o 



J"(3) + 0.85 J?(c) 



/'o 



(si; 



D" 



Zz/ + X. 



Determination of ikL 



The purpose of this section is to find out the value of the c-quark "hard pole mass r 
(see sect. 2), M c , using the experimental data for the semileptonic width |25|], [[26J1 : 



T 8l (D + ->lv + X) = (1.06 ±0.11) - 10" 13 GeV . (82) 

The expression for T s i can be represented in the form (see || and sect. 3a): 



r,z(£- 



1 si 



&n -*■', 



< D\cc\D> 



o 1-rad 



2M, 



D 






* °lept 



(83) 



1 si 



G\Ml 

192 7T 3 



l rad 



2 a s (M?) m 2 

3 7T /0l M2 J 



M) 



where 5\f vt in Eq.(83) is the contribution of SL^l into T s i, and the explicit form of 



f {x) can be found in p7| , p8fl . We use below Q: 

Q /l/f2 _ /Lf2 

m s ~150MW, <(/) c >~ 0.3 GW 2 , A G ~ - g * - g 

4 Mr 



55) 



The value of M c can be determined now from a comparison of Eqs.(82) and (83). 
Because the dependence of T s i on M c is highly nonlinear, it is more convenient to pro- 
ceed in an opposite way. Namely, let us show that Eq.(83) reproduces the experimental 
value at M c ~ 1.65 GeV. We have: 



aJM 2 ) ~ 0.310, ^ ~ 0.938, z t ~ 0.967. 



W) 



A G ~0.15, A^~ 0.055, / ~3.25. 



*7) 



5 From our viewpoint, the value: < p 2 >~ 0.6 GeV 2 available in the literature [£9| is overestimated. 
Let us recall (see Eq.(42)) that the mean value of the vacuum quark 4-momentum squared is: < 
— k 2 > = 4/3 < k 2 > ~ OAGeV 2 , and the quarks inside the pion have their momenta on the 
average somewhat less than those of the vacuum quarks |50| . Let us point out also that < p 2 > enters 
here to < D\cc\D > only and plays no essential role. 
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T° sl ~ 2.80 • 10" 13 GeV , 



I rad ~ 0.786 , - — ^r — - ~ 1.02 
2 M D 



1-2^A G 



0.691. (89) 



Substituting all this into Eq.(83) one obtains: 

lept 
(c) 



T sl (D + )^r° sl (0.52 + 5l:l t ). (90) 



The quantity 8\ e L is calculated in sect. 11 and is: 

jf(c) 

J Zept 



<fi~-0.13. (91) 



So: 

r si ( J D+) ~ 0.39 r° sl ~ 1.09 • 10~ 13 GeV , (92) 

in agreement with Eq.(82). 

5. Mass formulae. Determination of M b and A 

To find out the value of the b-quark "hard pole" mass, Mb, we use the mass formula 
which has the form (see, for instance, |17| and sect. 2) : 



M D = M c + A+^m 2 + o(^V (93) 



. 2 1 < D(p)\cAc\D(p) > 1 „ 3 /lrf „,2n 

6m S 2 2Af c 2 < p1 > "8 (M - - M °>' (M) 

A = {iDf - \g,a^Gl^. . (95) 

As it was pointed out in sect. 2, the difference between the "hard pole mass", M c , and 
the pole mass, M p , can be absorbed by a redefinition of A appearing in the matrix 
element: 

1 / A 2 \ 

< D\ H Ught \D>=k [l + 0{-&g-) , (96) 



2M D ' " y " M u \ y M 2 



or, equivalently, through the trace anomaly: 



1 . m0(«)^ lrw_A , 2 ^ , ^ A QCD^ 



< D\ %zr Gi v \D>=A + —— + O(^f) , (97) 



2M D ' 2na ^ v ' " M c M 2 
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and by a redefinition of the next term, (1/4Md) < D\ cA c \D >= 5m 2 — > 5m 2 , etc. 
Turning now to concrete numbers, we use in Eq.(93): 

Af c ~ 1.65 GeV, M D ~ 1.867 GeV, 
< f c >~ 0.3 GeV 2 , ( M 2 D , -Ml)~ 0.543 GeV 2 , (91 

and obtain: 



1.867 GeV ~ 1.65 Gel/ + A + 91 MeV - 123 MeV 



(99) 



A ~ 250 MeV 



(100) 



Using now this value of A and: 



' =* - 3 < Ml*-Mi)-~ O 



M ^ Mb + A+ 2M b <P » > -8M b ^ 



,A7 



M 2 



(101) 



M B -5.28 GeV, < f b >~ 0.3 GeV 2 , (M|* 
one obtains from Eq.(lOl) Q 

M 6 -5.04 Gel/. 



M|) -0.488 GeV 2 , (102) 



(103) 



6 . B —* I v + X . Determination of I V C J . 



The semileptonic decay width " B — > e v + X" is obtained from Eq.(83) by evident 
replacements. We have: 



J- =; 



G^M^I 2 

192 7T 3 



~ 1.19- 10~ 13 GeV 



V cb ' 



0.040 



(104) 



a s (M 2 ) -0.204, 2 D ~ 0.460, z x ~ 0.635, / ~ 2.46 , (105) 



< B\bb\B > 



2M 



1.00. 



B 



1-2^A G 



0.96 



(106) 



'■rad 



2 a s (M 2 ) M 2 C , 
3 7T M M 2; 



0.893 



D Really, it is clear that we need not A at all to determine Mb . 



(107) 
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So, (see Eq.(83): 



T(B ^ ev + X) ~r° sl (0.394 + 5® t j . 



The quantity 5\J t is calculated in sect. 11: 



8^1 ~ - 2 • 10- 3 , 
and is negligble. Therefore: 

r(B ->■ e z/ + X) ~ 0.47 ■ 1(T 13 GeV 
If we take, for instance, |2(J [] : 

-12 



Kfc 2 



0.04 



r(B) = (1.6 ± 0.04) • 10~ 12 s , r tot ( B ) = (4.1 ± 0.1) • 10~ 13 GeV 



'108) 



(109) 



(110) 



:m: 



Br(B ^lu + X) = (11.4 ± 0.5) 



then: 



(112) 



r (B -»• / 1/ + X) = (0.47 ± 0.03) • 1(T 13 GeV 
and comparing with Eq.(llO) we obtain: 



V cb \ ~ 0.040 



Br(B ^ lu + X) 



Br(B^lv + X) 

11.4% 



11 x 



nV2r 1.6-10- 12 s ll/2 



1.6- 10-^s 
r(B) 



i;i3) 



r{B) 



(114) 



7. Calculation of fa , /b 

There is a large number of papers dealing with the calculation of the decay constants 
fjj and Jb with the help of the QCD sum rules. We point out here only early papers 

0,0 



|3"3"|| , and the paper [^3| close in spirit to this work, were it was proposed to 
use the correlators of chiral currents in the QCD sum calculations. 

7 We use the LEP data for the b-quark lifetime and Br(b — ► eu + X). When comparing these with 
the T(4S) data, one sees that the absolute value of T(B — * ev + X) ~ T(b — > e^ + X) is the same, 
while the LEP data give a smaller value for the total decay width. We prefer to use the LEP value 
for the following reasons. It gives the weighted average of the b-hadron widths. Because (see below) 
the B-meson widths are practically the same, and there is all the reason to expect that the b-baryon 
total widths are somewhat larger than those of the B-mesons (with the semileptonic widths being the 
same), the admixture of the b-baryons can only increase the weighted total decay width. So, the LEP 
data can be considered as giving the upper limit for the B-meson decay width and the lower limit for 
the semileptonic branching. 
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Due to reasons which are explained below in detail (see also J34|), our approach here 
rests heavily on the use of the chiral currents correlators. One of the main difficulties 
which prevents the calculation of reliable results from the correlator of the pseudoscalar 
currents, is the very large radiative correction to the Born approximation for this 
correlator [35], [36], [3?]. To avoide this difficulty, let us consider the following correlator 
of the chiral currents: 



KM 2 



/ / dx e iqx < 0|T {Q(x)i(l + 75) q(x) , g(0)z( 1 + 75) Q(0)} |0 

fdxe iqx < 0\T{P{x)P + {0) - S{x) S + {0) } |0 > , 
P{x) = Qi^ 5 q(x) , S(x) = Q(x) q(x) . 



>z 



(115) 



Because Ki(q 2 ) is the difference of the pseudoscalar and scalar current correlators, 
the pure perturbative contributions cancel completely in all orders of the perturbation 
theory (in the chiral limit). On the other hand, there appear additional scalar state 
contributions in the spectral density. Let us emphasize, however, that the mass differ- 
ences between the lowest lying pseudoscalar and scalar resonances are sufficiently large 
nevertheless, both for the D and B mesons. So (after the Borel transformation), the 
scalar meson contributions are sufficiently suppressed in the sum rules. 

The diagrams giving the main contribution to Ki(q 2 ) are shown in fig. 5. The 
spectral density (for the D-meson) has the form: 



SK(s)=r 2 D 5(S-M 2 D ) + 



r D = f D M 2 D /M c 



(116) 



The sum rules obtained in a standard way have the form 0: 



r 2 D {M 2 ) = 2M C < 0\qq\0 >^ $ C (M 2 ) , 



(117) 



$ C (M 2 



exp 



M 2 D 



Ml 



M 2 



m 2 M 2 ( i ^M 2 ' 
4M 4 \ M 2 



:ii8i 



m„ 



< 0|g|a^g^-gjO >jp _ < 0|g(?£> M ) 2 g|0 >^ 



< 0|gg|0 > 



M 2 



< 0\qq |0 > 



M 2 



;ii9) 



In Eqs.(117) and (119) M 2 ~ lGeV 2 (see below) is the normalization point of all 
operators and of the D-meson residue, Td , and the number m 2 is also determined at 
this point. We will use for it the standard value determined previously |23| at this 
scale: m 2 ~ 0.8 GeV 2 . 

{Fig. 5} 



8 We neglected in Eq.(118) the small contribution ~ < qq > 2 . The anomalous dimensions of the 
operators q~D 2 q and qq are respectively: (— 2/36 ) and (4/6 ) . 



Figs. 5a, 5b, 5c The diagrams contributing to the sum rule Eq.(118) 

The power corrections due to the figs. 5b, c diagrams do not exceed ~ 35% of the 
fig.5a contribution in the region 0.8 GeV 2 < M 2 < 1.5 GeV 2 and, on the other hand, 
the values of M 2 are sufficiently small for the contributions of higher states in the 
spectral density to be exponentially suppressed. The function $ C (M 2 ) varies only 
slightly in this region and the characteristic value of M 2 (at the extremum of $ C (M 2 ) ) 
is: 

M 2 ~ 1.15 GeV 2 , $ C (M2)~1.8. (120) 

Because we need in what follows the value of fo at the low normalization point /z , we 
renormalize now both r 2 D and < qq > to this point and (because they have the same 
anomalous dimension) obtain: 

r 2 D (^ 2 )~2M c <0\qq\0>, o 1.8 . (121) 

Finally, to obtain the concrete answer we use 

<0\qq\0>^ o ~ (0.25 GeV) 3 , \x 2 ~ 0.5 GeV 2 , (122) 

for the value of the quark condensate at the low normalization point. So, 

r D (jj, 2 o = 0.5 GeV 2 ) ~ 0.30 GeV 2 , f D (/j 2 ) ~ 144 MeV , 

fo(M 2 ) = (^j0y) fn(nl) * 165 MeV . (123) 

For the B-meson, the corresponding sum rule has the form: 

r 2 B (W) =2M b < 0\qq\0 > m $ 6 (M 2 ) , (124) 



, n N {Ml -Ml 

$ fe (M 2 ) = exp I h '' 



1 r< M Ul 2 M ^ 

4M4 M 2 



M 2 
The factor k in Eq.(125) (M 2 ~ AGeV 2 , a s (AGeV 2 ) ~ 0.284 



(125) 



/ a s (W)\ U/25 ( a s (M 2 ) \ 14/27 , x 

is due to different anomalous dimensions of the operators qig s aGq and qq . 

The corresponding interval (see above) of M 2 is: 3 GeV 2 < M 2 < 5 GeV 2 , the 
function $^(M 2 ) varies only slightly in this interval, and the characteristic value of 
M 2 (at the extremum) is: 

W^AGeV 2 , $ b (M2)~1.53. (127) 

19 



Proceeding now in the same way as for the D-meson, we obtain: 
r B (fi ) ~ 0.49 GeV 2 , f B (jjt ) ~ 89 MeV , 

MiU) = GS)"* t^Sf'^ ^ 113MeX/ - (128) 

It is of interest to compare the above results with those obtained in the static 
limit: M q —* oo. Let us define the constant f (a static analog of fjj) as: 

M/io) = -/^/o, (129) 



M D ~M C + A, M 2 = 2M C E, A ~ 250 MeV. (130) 

Then, in the limit M c — > oo, the sum rule Eq.(118) takes the form |: 

f o ^2<0\qq\0> llo( j ) (E) 1 <f>(E) = e K ' E (l - ^U . (131) 

The corresponding interval (see above) of E is: 400 MeV < E < 700 MeV, the 
function <p(E) varies only slightly in this region, and the characteristic value of E (at 
the extremum) is: 

E~ 550 MeV, </>(E ~ 550 MeV) ~ 1.32 . (132) 

Substituting Eq.(132) into Eq.(131) we have: 



/ o ~0.20GeV 3/2 , f D {M c 



Os(jJ%) 

a,(M c 2 



2/9 



/ ~ 180 MeV. (133) 



/M 



c 



Comparing Eq.(133) with Eq.(123) we see that the result obtained for f D in the static 
limit does not differ considerably from those obtained for the real value of the c-quark 
mass. 

Although we succeeded in using the correlator Eq.(115) in which both the Born 
contribution and all radiative corrections to it are absent, for the calculation of the 
decay constants fjj and fs, the analogous trick, unfortunately, turns out to be useless 
for the calculation of nonfactorizable contributions to the matrix elements of the 4- 
quark operators which we will need in what follows. The reason is as follows. In 
the sum rules for the correlator Eq.(115) the leading contribution is proportional to 
the known quark condensate < qq >, and the main correction is proportional to the 
condensate < qaGq > which is also known. However, if we will try to calculate the 
matrix elements like (see below): < D\cT ^q ■ qT u ^c\D > in an analogous way, 
then even the leading contribution is expressed through the high dimension condensate 



Really, we drop out power corrections but keep the value of rn^ fixed. 

20 



< qqqqGG > which is unknown. Therefore, the only way to calculate such matrix 
elements with the help of the QCD sum rules is to use such interpolating currents that 
the leading contribution is pure perturbative, and the leading corrections are expressed 
then through the known low dimension condensates like < G 2 > , etc. 

As it was pointed out above, to deal with such sum rules we have to overcome 
in some way the difficulties originating from very large radiative corrections. We will 
describe now the way we used here to deal with this problem and which is used then 
in sect. 9. 

Let us return to a consideration of the decay constant fn and consider now the 
correlator: 

K 2 (q 2 )=ijdxe^<0\T{Q(x)(l + lb )q(x), g(0)(l - 75 )Q(0)} |0 >= 

% J ' dxe iqx < 0\T [P{x) P+(0) + S{x) S + (0)} \0 > , (134) 

The sum rule obtained from Eq.(134) in a standard way has the form (in the chiral 
limit): 



2/ , 3 fSodS(S-MiY M 2 D -S 

r D<J*o) ^73/ 5 exp ' 



'* '"> - r 1 + i-M).,(M) + 0(aM , 5) 



X(z) 



4vr 2 Jm S * M 2 



9 3 z 1 z 

- + 2 l(z) + \og{z) iog(i -z) + - log J-— + log yzt z ~ z log T^z~ + 



z , 1 3 . Ml 

i log - - 7 log — 

1 — z z 4 jj, a 



o 



L(z) = - r^log(l-t). (136) 

jo t 



Let us emphasize that the chirality odd condensates like: < qq > , < qaGq > etc., give 
no contribution to the correlator Eq.(134). Q 

Let us try now, as a first approximation, to neglect all radiative corrections in 
Eq.(135) and let us choose the value of S (the effective parameter which models the 
beginning of the perturbative continuum in a given correlator) in Eq.(135) to obtain a 
fit in M 2 . There is a good fit in the standard region (see Eq.(120) ): 0.8 GeV 2 < M 2 < 
1.5 GeV 2 at S = 3.81 GeV 2 , which gives: 

r D (p ) ~ 0.096 GeV 2 , f D (ji ) ~ 45 MeV . (137) 

For the B-meson, there is a good fit in the standard region 3 GeV 2 < M 2 < 5 GeV 2 at 
S = 28.8 GeV 2 which gives: 

r B (fx ) -0.19 GeV 2 , f B (fi ) ~ 35MeV . (138) 

It is seen that the Born approximation to the sum rule Eq.(135) gives (at M c = 
1.65 GeV, M b = 5.04 GeV) very small values for f D , f B (compare with Eqs.(123), (128)). 

10 We neglect the small contribution to Eq.(135) due to < G 2 > and < qq > 2 . 

21 



Let us account now for the one loop correction to the Born approximation in the sum 
rule Eq.(135) and make anew the fits in M 2 . We will obtain good fits both for D- and 
B-mesons in the standard regions of M 2 - values at S = 3.85 GeV 2 and S = 28.9 GeV 2 
respectively with the results: 

r D (fi ) ~ 0.149 GeV 2 , f D (pt ) ~ 71 MeV . (139) 



r B (fi ) ~ 0.324 GeV 2 , f B (n ) - 59 MeV . (140) 

It is seen from a comparison of Eqs.(139),(140) and Eqs.(137),(138) that if, in the Born 
approximation, we obtain the values of fo, fs which are ~ 3 times smaller than the 
right values Eqs.(123), (128), on account of the first radiative correction the results 
increase ~ 1.6 — 1.7 times. It is clear that, in such a situation, one needs either to 
account for all radiative corrections, or to use some trick in the hope that it can help 
to account effectively for a summary effect of all radiative corrections with reasonable 
accuracy. We have chosen the second way, of course. 

The large radiative corrections in the correlator Eq.(134) look reasonable as they 
tend to increase the Born contribution which is too small. On the other hand, these 
large corrections indicate that the values of the parameters entering the spectral density 
in Eq.(134) are not chosen properly. In the given case, the only parameter is the quark 
mass, M c , (we used everywhere above the "hard pole masses"), and the correlator 
Eq.(134) is very sensitive to the precise value of the quark mass (with the D- meson 
mass fixed). The right hand side of Eq.(135) increases greatly when the quark mass 
decreases. 

Therefore, it looks natural to try to describe the main effect of radiative corrections 
by using the effective quark masses: /i c < M c , fib < Mb. So, let us define: 

M c = fijl + ^4^- C D + 0(a 2 )) (141) 



(and analogously for Mb), and let us express M c through fi c at the right hand side of 
Eq.(135). We obtain: 



... , ., f SodS(S-fi 2 c ) 2 \M 2 D -S' 



4ir 2 Jul S H M 2 



l + l^^X(f) + 0(^)(lJ2) 
on b 



X{z)=X{z)-\c D -^- z . (143) 

Let us try now to find a value of ji c which will give the answer for fu close to the right 
one, Eq.(123), and for which the radiative corrections will remain reasonably small at 
the same time. It is clear that there is no guaranty that it is possible to succeed in this 
way. 
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Nevertheless, the results are very encouraging. For instance, let us suppose that 
the radiative corrections are reasonably small (for a properly chosen value of /j, c ), put 
X(z) = in Eq.(142) and find the value of /j, c which reproduces the right answer, 
Eq.(123). For: jj, c = 1.40 GeV, S = AAGeV 2 one obtains a good fit in the standard 
region 0.8 GeV 2 < M 2 < 1.5 GeV 2 with the result r D (fi ) — 0.3 GeV 2 (compare with 
Eq.(123)). To elucidate the role of radiative corrections in the sum rule Eq.(142), let us 
calculate it now in the same region and with the same papameters: \x c = 1.40 GeV, S = 
4AGeV 2 , and Co = 3.87 (see Eq.(141)). The fit is not optimal but sufficiently 
good, and when the right hand side of Eq.(142) is taken at the characteristic value 
M 2 = 1.15 GeV 2 it gives: r D (fi ) — 0.29 GeV 2 , in good agreement with Eq.(123). 

For the B-meson, the situation is analogous but slightly worse. Namely, let us 
choose /j,b = 4.84GeV and put X(z) = 0. We obtain then a good fit in the standard 
region: 3 GeV 2 < M 2 < 5 GeV 2 at S = 30.0 GeV 2 with the result: f B (fi ) - 89 MeV, 
which reproduces Eq.(128). Accounting now for the radiative correction, using \Xb = 
4.84 GeV, S = 30 GeV 2 , and C B = 0.895 (see Eq.(141)), and taking the right hand 
side at the characteristic value M 2 = 4 GeV 2 , we obtain the result for tb{^o) which 
is ~ 15% higher than the right value, Eq.(128). Fortunately, because the power 
corrections to the B-meson decay widths are small (see below), such accuracy will be 
sufficient for our purposes. 

In summary, it is possible to reproduce the right values of fu and fg by neglecting 
the radiative corrections in the sum rule Eq.(142) and using the effective quark masses: 

H c = lA0GeV, i2 b = 4.84GeV. (144) 

The residual effect of radiative corrections remains reasonably small in this case. 

As for the correlator Eq.(134) and the decay constants fo and /#, there was no 
need to perform all the above manipulations because the answer was obtained previ- 
ously in this section. As it was pointed out above, our real purpose is to calculate 
more complicated matrix elements of the 4-quark operators (see sect. 9). And our main 
assumption is that for the correlators used below in sect. 9 it is possible to obtain rea- 
sonable results by neglecting radiative corrections and using \x c and \x\, instead of M c 
and Mb in the corresponding spectral densities. 

8 . Difficulties with naive estimates. 

Let us turn now to T° nl defined in Eq.(25). All the quantities entering are known 
(because the radiative corrections are known here for m s = only and we neglect below 
the SU(3)-symmetry breaking corrections, we put m s = in nonleptonic calculations): 

for 2 4- G 2 \ ( G 2 - G 2 \ 

[ % - j^l-43, (1-2 As) -0.7, 4 ^ - ^ j A G ~ 0.36 , (145) 

<D }f} D> c^l.02, I rad ^l.05, T o = ^$^8A.l0^GeV, (146) 
I Md 64 7T" 3 
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T° nl ~ r o [1.43] md [1.02U [1.05] md ^ 1 - [0.3] ctG j + [0.36]^ j ~ 1.59 r o . (147) 

Let us try now to obtain a rough estimate of the matrix element < D + \SL e ff\D + >, 
see Eq.(45), by putting: A ~ Pp in L s , L u (see fig. 3, Pp is the D-meson momentum), 
A 2 ~ P c 2 ~ M 2 in L d and using the factorization approximation [0], |3"8| . Then, L u and 
L s give zero contributions and: 

< D+(p)| cT^ d ■ dT v c \D + (p) > Mo ~ PflPu f 2 D (/i ) , (148) 

< D+| L d |£>+ >£*"-~ -1.1 M 2 / 2 (^ Ml , (149) 

2 / 2 (/i )M D " 



^r/actor(-D 



-1.1- 16 7T' 



r o ~-i.50r o . (150) 



M c 3 
The contribution from Lp^y in Eq.(51) is obtained by using: 

< D\ j; • c 7p (l ± 75 ) y c \D >£**- -- rU»o) { 1 - ^ J , 



(151) 



7T 



^(l-.Wl-^U-O.lS. (152) 



r o ~ 9 M c 5 "V ,co /^ 2M| )/ 

On the whole, one obtains: 

Y nl (D + ) ~ (1.59 - 1.50 - 0.15) Y ~ -0.06 r o , (153) 

which does not make much sense. It is clear that the above approximations are too 
rough and some estimates are essentially wrong. It is the purpose of subsequent sec- 
tions to improve the above described naive estimates. 

9 . Nonfactorizable contributions: gluon condensates. 

We calculate in this section the nonfactorizable gluon contributions to the matrix 
elements of the 4-quark operators. Let us begin with the operator 0^ u = cT^(X a /2)q ■ 
qT u (X a /2)c, which gives zero matrix elements in the factorization approximation. We 
want to calculate the contribution of the fig. 6 diagrams with the help of the QCD sum 
rules. 

{Fig. 6a} 

Fig. 6a The nonfactorizable gluon contribution to the weak annihilation 

{Fig. 6b} 
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Fig. 6b The nonfactorizable gluon contribution to the cross weak annihilation 

{Fig. 7} 

Figs. 7a, 7b The diagrams for the meson transition into a current with an emission 

of a nonperturbative gluon 

Practically, however, it is more convenient to calculate the nonperturbative gluon 
emission amplitude, fig. 7, and to obtain then the contribution to the decay width by 
averaging its modulus squared over the gluon field fluctuations in the QCD vacuum. 

To calculate the fig. 7 amplitude, we replace the D-meson by the interpolating 
current and consider the correlator: 

T£ = i fdxexp{ipx} < lgl \T J°(x) J P (0) |0 > Mo , (154) 

j;(x) = q (x) r M y c{x) , J P (0) = c (0)i( 1 + 75 ) 9(0) . (155) 

Calculating the contributions of figs. 7 a, b diagrams to the discontinuity of Tf in p 2 , 
integrating it from the threshold up to S with the weight exp{— p 2 /M 2 } and equat- 
ing to the D-meson contribution, we obtain the sum rule for the D-meson transition 
amplitude into the current J" with an emission of a nonperturbative gluon 0: 

N » - (2 ^)vX) {* dU lG + G% ^ Pa ' (156) 



T; 



;f -p{%^} ; 'o=L -"■„ " c ' ^^ 7 ^^. (157) 




"/if „ tfiv Act "Act • 

(In accordance with the discussion in sect. 7, we replaced M c by /i c in the spectral 
density, see Eq.(144), supposing this accounts for the major effect of radiative correc- 
tions). 

Unfortunately, the attempt fails to obtain a good fit for Iq in the standard region 
of M 2 : 0.8 GeV 2 < M 2 < 1.5 GeV 2 , and with So in the reasonable vicinity of its 
optimal value S = AAGeV 2 (see sect. 7), because Iq is nearly independent of S . For 
instanse, I G varies as: 2.0 > I G > 1.2 in this region at S = AAGeV 2 . Therefore, the 

n The term Iq in Eq.(156) is due to the fig. 7b diagram and is parametrically smaller than those 
from fig. 7a. 
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best we can do is to take for I G its value at the characteristic point M 2 = 1.15 GeV 2 
(see sect. 7) : 



T G ~ I G (M 2 = 1.15 GeV 2 ) ~ 1.4. 



(158) 



Analogously, 0.4 > I G > 0.3 in the standard interval of M 2 and at 5 = 4.4 GeV 2 . So, 
we take: 

I G ~I G (M 2 = 1.15GeV 2 ) -0.33. (159) 

On the whole, we estimate: 



N ~ 
J V — 



4vrM c r7 s 



27r) 4 r D (/v) 



{l.4zG» a +0.33G« a } Pa . 



(160) 



Taking now the product N^ iV+ and averaging the gluon fields over the vacuum, we 
have finally Q: 



<D(p)\ 



Y 



cT ^Y q ' qTu ~ c 



/to 



\D(p)> 9 ^[^-gJc G , 



(161) 



MlM 2 



a 



48 7r 2 rf,(/i ) 



<0|— Gf,„|0> 



7T 



^ti^ I 



/g) 2 



/ G ) 2 



0.48 ■ 10~ 2 GeV 4 . (162) 



The situation is qualitatively the same for the B-meson (I G varies as: 0.46 > I G > 0.35 
in the standard region 3 GeV 2 < M 2 < 5 GeV 2 , S = 30 GeV 2 ), and proceeding in 
the same way we obtain I G ~ 0.4, I G ~ 0.03 , so that: 



•tb 

a 



C^~1.15-10-"Ge^ 4 . 



(163) 



Let us point out that the parametric behaviour of Iq is: Iq = 0(/i o /M c ) , Iq = 
0(/z 2 /M 2 ). So, the matrix element Eq.(161) which describes the non-perturbative con- 
tribution to WA (weak annihilation, fig. 6a ) and CWA (cross weak annihiltion, fig. 6b 
) behaves as: 0(/^M c ), i.e. in the same way as the factorizable contribution Eq.(148). 
This gives a relative correction to the decay width: 5T/T = 0(^/M|), as it should be. 

In order to understand to what extent the factorization approximation is good, let 
us compare Eqs.(161)-(163) with the factorizable matrix elements: 

< D{p)\cT, q .qT v c\D(p) >£*- {^1 fUn )Ml ~ (^) 7.2 • 10~ 2 GeVfai) 



<B{p)\bT,q-qT v b\B{p)>^ 



P^Pv 
p 2 



f B (»o)M_ 



B 



P^Pv 
p 2 



22.0 • 10" 2 GeV 4 



12 



It is implied, see figs. 6, 7, that the quark flavour in Eq.(161) is such that it is a valent one. 
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It is seen that the nonfactorizable contributions are ~ 15 — 20 times smaller than the 
corresponding factorizable one, so that the factorization approximation works very 
well even for the D-mesons. 

It seems therefore that there are no chances to change essentially the results ob- 
tained in the previous section. As will be shown below, this is not the case really for 
the following reasons: 

1) Although the matrix element Eq.(161) is small, it enters SL e ff with much larger 
coefficients in comparison with the factorizable contributions; 

2) As will be shown in sect. 12, the characteristic value of A 2 in La, Eq.(31) is: 
< A 2 >~ 0.35 M^, while those of A 2 in L u , L s is: < A 2 >~ M 2 D , and this effect 
suppresses strongly the large factorizable contribution, Eq.(164). 

Let us perform now an estimate of the nonfactorizable contribution to the matrix 
element of the operator S iliV = cT^q ■ qT v c. For this, let us consider the correlator: 



U 



j dxexp{ipx} < 2gl\TJ IJi (x)Jp(0) |0 > flo , 



(165) 



Mx) = q{x)T fl c{x), J P (0) =c (0)z(l + 7s ) g(0) 
and calculate the fig. 8a contribution. We have: 
8iM c g 2 f dk 



n 



\k k v g^p kp k v k^ 



2<kY J k 8 [(p + k) 2 -M 2 ] 
One obtains from Eq.(167) with logarithmic accuracy: 

2 . r . ///:' \ d 1 



2?~fi M -M c log 



r'max 
, r'min , 



dp 2 p 2 - M 2 



pX ^\u 



(166) 



(167) 



(168) 



a 



n 



(jfiX&pXPp ~7Pn( jr p\G p \ 



(169) 



Proceeding now in the standard way, one obtains the sum rule for the D-meson tran- 
sition amplitude into the current Jp and two gluons: 

JM 2 D -MP 



M^~ 



2 1 
Q„ ~M C —— 

3 r D (fi c 



r'max 



log ^ — exp 



, r"min 



M 2 



M 2 



+ ■■ 



(170) 



As we need here a rough estimate only, let us put: M 2 = lGeV 2 in Eq.(170) (see 
sect. 7). Besides, because both \j? ma3 . and /i^ in in Eq.(170) remain finite at M c — > oo, 
we put: 



log 



r'max 
. rmin , 



1. 



(171) 



13 Because we can expect beforehand (see the above discussion) that the nonfactorizable correction 
will be small, we confine ourselves to the main contribution from the fig. 8a diagram, neglecting even 
smaller contributions from those diagrams where gluons are emitted by the c-quarks. 
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Therefore, 



M M ~ (7.9GeV- 3 ) Q„. (172) 



Now, taking the product M M Mj~ and averaging the gluon fields , one has: 

-6 



<D(p)\cT ll q-qT l/ c\D(p)> nonfact ^[Q2GeV' b ) < 0|n„n„ |0 > . (173) 

Unfortunately, this 4-gluon vacuum condensate is unknown. So, let us obtain first an 
estimate in the factorization approximation, see fig.8b. 

{Figs. 8a, 8b} 

Fig. 8 The nonfactorizable contribution to the matrix element < D\S^ U \D > 

One has: 

<0\n ix n v \0>~-±g lil ,M* ) <0|^Gjj0> 2 +O(p M ^), (174) 

so that: 

< D(p)\ cT^q ■ qV u c\D(p) > nonfact ~ g„ v (0.5 • 10~ 4 GeV 4 ) , (175) 

(and the term ~ p^p v is of the same order). Comparing Eq.(175) with Eqs.(164) and 
(161) we see that the nonfactorizable correction is very small in this matrix element. 
Therefore, even if the approximations made above give the right order of magnitude 
only, we can safely neglect this correction. 

10 . Nonfactorizable contributions: quark condensates. 

The contribution < D\5Lpnv\D > wa s calculated in sect. 8 in the factorization 
approximation. Because it is not large by itself and nonfactorizable corrections are 
also small (see sect. 9), it does not make much sense to account for them in this matrix 
element. So, we will use for it the expression Eq.(152). 

There are, however, analogous corrections of the "penguin" type in the matrix 
elements of the operators L u , Lj, L s in Eqs.(45)-(48), see figs. 9, 10, and we proceed 
now to their calculation. 

{Fig. 9} 

Fig. 9 The nonvalence nonfactorizable penguin contribution (NV) 

{Fig. 10} 

Fig. 10 The valence nonfactorizable penguin contribution 



For a calculation of the diagrams in figs. 9, 10 (see also fig.4) let us consider the 
correlator: 



V = % \ dxe^ x i / dye-*™ < 0|T J+(s) S M „(0) J F (y)|0 >, 



;i76) 



Jp = ci(l + 7 5 )g, Jp = gi(-l + 7 5 )c, S^ = cT^ • ^T v c. 



177) 



where q and ■?/> are light quark fields. 

One can neglect the co-ordinate dependence of the light quark fields in Eq.(176) for 
calculations with logarithmic accuracy, and obtain: 



T 



2M r 



2 M r 



v\-Ml \pt-Ml 



K u „ + - 



(178) 



Kpv =< 0|g(0)r>(0) • $(0)r„g(0)|0 > . (179) 

Proceeding in the usual way one obtains now from Eq.(178) the sum rule: 



< D(p)\S^(0)\D(p) > pen9 ~ K t 



fl V 



2M C 



D c 

— 72 — 



e "i + ... 



2M C 



■e M 2 +... 



,(180) 



where the dots denote power corrections in 1/Mf , 1/Mf . 

Let us consider now the contribution of various penguin- diagrams to the matrix 
element < D\5L e ff\D >, (see Eq.(45)). In order to check the normalization, let 
us start with the contribution of fig.4. Separating in K fll/ , Eq.(179), the penguin 
contribution from the operator if) if): 



i/ji/j 



l-Vco 2/b ° f Xa 



ip J a P m , 



peng Q I 2 

and factorizing in the standard way the 4-quark condensate, we obtain: 

< D\S, V \D > ftgA ^ A giiv (i _ ,,-2/fc-) 



:isr 



2M c <qq> 
rD 



e M i +... 



2M c <qq> ^ 



td 



e M 5 + .. 



:i82i 



On the other hand, to calculate this contribution there is no need to use the sum rules 
at all because, after the replacement Eq.(181), this contribution to the matrix element 
can be obtained directly through factorization: 



' i „— 2/b \ , 

< d\ s, p - ( — ^_ ) cr,- 7p r> • j; \d > figA ^ - <w (1 - r 2/6 °) ^(183) 



6 



~2 
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Comparing with Eq.(182) we see that we simply reproduced the sum rule used before 
(see Eqs.(117),(118)): 



2M C < 0|gg|0 > e" 



{2 D - M c 

M' 2 



[1 + 



D 



Therefore, we can rewrite Eq.(180) in the form: 



<D\S liV \D>>. 



D 



< 0\qq\0 >* 



K, 



fll/ ■ 



;i84) 



(185) 



It is not difficult to obtain now other penguin contributions from Eq.(185). Namely, 
the contribution of fig. 9 is obtained by separating out the penguin contribution from 
the operator [qq] in K^ v : 



qq 



Peng 



a s\^o) , { Umax 

2vr g 



r^min , 



$)*• 



;i86) 



which gives after the standard 4-quark condensate factorization: 

2 



< D \ [<VU \ D >fig.v- 2j9t*vr D (jio) 



2-n 8 



H"min , 



A, 



(187) 



where A s = 1 if the ip in S^ v is u- or d-quark, and A s = (< ss > / < uu >) 2 ~ 0.64 
if ip is the s-quark. Let us point out that, unlike the standard penguin contribution, 
fig. 4, which contains the factor: (1 — i]~ 2 / bo ) ~ (a s (fil)/2ir) log(M c 2 //i 2 )), the fig. 9 
contribution contains as an upper cut off the quantity \J? max which remais finite at 
M c -> oo (see also Eq.(171)).Q 

If we consider now the operator O IMU = [c^-T^ip -ip^-T^c] instead of S^ u , Eq.(177), 
then we have clearly: 



<D\O liV \D> 



fig-9~~ 



~ 7 < D\Sfj_ u \D >fig.9 



For the operator M ^, however, there is an additional (valence) penguin contribution 
originating from the diagram fig. 10 (plus the mirror one). It can be obtained easily 
from the relation analogous to Eq.(185): 



<D\0^ u \D>: 



D 



< 0|gg|0 > 



p 



(189) 



A" 



P^=<^\qT^'—q-qV v '—q\Q>, 



A^ 
Y 



(190) 



14 It is clear that the contributions of figs. 4, 9 are nonvalence, i.e. the same for all Z? 0,± ' s -mesons; 
they only shift the position of the " decay width centre" F° z . 
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by separating out the penguin contribution: 



qT,-q 



peng 



a s{Ho) i / Umax 



67T 



, AV, 



ja 



(191) 



and factorizing the 4-quark condensate. We obtain (the factor 2 accounts for the mirror 
diagram) : 



< D \ [CWU \D >fi 9 .w^ ^9^r 2 D {fi ) 



a s\^o) i / fJ'max 



2vr 



. Hrnir. 



(192) 



Let us point out finally that, within the logarithmic approximation, we can put: A ~ 
A ~ P c ~ M c for the contribution of fig. 4, and A ~ A ~ Pjj for those of figs. 9, 10 . 

For the B-mesons, it is sufficient to replace: r£>(p ) — > rsiPo) m the above formu- 
lae Eqs.(187), (188), (192) to obtain contributions to the matrix elements of L u , L^, L s 
in Eq.(61) (for the B ±,0,s mesons there are no contributions from L c , L c , L c in Eq.(61)). 

11. Corrections to semileptonic widths. 

We are ready now to explain how the value of 5 lept (D + ) used in sect. 3 was obtained. 
Let us recall (see sects. 8, 10 and fig. 9) that: 



2 / < ss > 

< D + (j>)\ cr M s ■ s T u c \D + (j>) > IMo ~ —g^ Po — — — 

27 V < uu > 



rl(Vo) , 



A<\ 



< D + (p)\c-r^s-s-r u c\D+(p) > Mo ~ -- < D + {p)\cV^s-sV v c\D + {p) > Mo , 



Po 



y 

2vr 



1 
6 



log 



r^max 
v r"tnin , 



0.1 



< 0|ss|0 > 

< Ol-u-ulO > 



0.8 



(193) 



< D+(p)\ cT^d-dT u c \D + (p) > Mo ~ p^p u f D {p ) + — g„ v Po r 2 D (p ) . (194) 
Therefore (see Eq.(55) in sect. 2): 



AF lept (D- 



2M D 



<D + (p)\AL l ^(p )\D + (p)> IIo : 



G 2 

"''' M 2 D p T co r 2 D {p )m o , ^ 



54ttM d 



< ss > z 

< uu > 5 



IK, 



IK 



c(/ 



0.656, (195) 



Ar ie P tp- 

■pfept 
1 o 



32 2 /£(/x )M 



-IT 



D 



Ml 



p T C0 ^ ~ -0.045. 



(196) 
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In addition (see sects. 3a, 8 and fig. 4) : 



r-plept 

01 PNV 



■ lept 



2M 



<D\L^ v ^ )\D>, o - 



D 



Lfp 



12ttM d 



^) Tm M?<D\ 



I CO lv± c 



r<2 



12ttM d 



: 2/9 )r™Mf 



^(3)-Cy 7/ ,(H-^75)c 

4 2/ J, M 2 \ 



\D > ~ 



(197) 



^v 64_ 2 /B(/io)M3 



-piepi 
1 o 



On the whole: 



-7T 



M« 



^l-^o 2/9 )r C o 1 



Ml 
2M D , 



-0.085. 



1 <D+\SL l e f f \D+> 



1 si 



2M 



D 



(198) 



-0.085-0.045)^-0.13. (199) 



For the B ±,0,s mesons (see sect. 3b), the term AL lept (/j, ) gives no contribution, and 
the analog of eq.(198) looks as: 



51%v 16 2 fjMM% 

3 



-plept 
1 o 



Ml 



-8/25 n, 
Vbc )TboV. 



2/9 (, Mj 

co I 2 M 2 



-2 • 10" 3 . (200) 



B, 



All the above described contributions are nonvalence. There is also a sizeable 
valence contribution to the D s leptonic width (see Eq.(56) and sect. 9): 



AT lept (D s 



G 



A-kMd 



F \y 12 



-0.245) T^ < D s \cT^— s ■ sT v — c\D s > Mo : 



G] 



AttM d 



\V CS \ 2 (-0.245) M 2 D {C C G -P, 



v ) > 



(201) 



^W|K,I>( 

1 o 



- 0M5 ^^o 



Py)^ 



-7%. 



(202) 
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A and A 



As it was indicated before, the values of A and A (in terms of the c-quark and 
spectator quark momenta) are clear from each diagram. Namely (P c is the c-quark 
4-momuntum, kj are the spectator quark momenta, Pd is the D meson momentum, 
Pd = P c i + ki — P c i + k 2 for the initial and final D mesons): 

a) A ~ (P c + ki) = Pd for the figs. 3a, 3b, 6a contributions; 

b) A ~ A ~ Pd (within logarithmic accuracy) for the figs. 9, 10 contributions; 
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c) A ~ A ~ P c ~ M c (within logarithmic accuracy) for the fig.4 contribution; 

d) A ~ ((-Pd — k\ — k 2 ) for the figs. 3c, 6b contributions. 

The differences between all of the cases above disappear in the formal limit Mq — > 
oo but, as will be shown below, they are of great importance for D mesons and are 
sizeable even for B mesons. 

Because the D meson momentum, Pd, is not an operator but a fixed number, we 
have to deal practically with the case " d" only. So, let us consider in detail the matrix 
element (see fig. 3c and sect. 3a): 

I D =< D(P D )\\ 2 cT^d ■ dT^c\D(P D ) >=<\ 2 D >< D\ cT^d ■ dT^c \D > , (203) 

where A = (Pd — k\ — k 2 ) and k\, k 2 are understood as the 4-momentum operators of 
the initial and final spectator quarks. 

As it was shown above (see sects. 9, 10), the factorization approximation works very 
well even for the D mesons, and the non-factorizable contributions are much smaller 
than factorizable ones. So, we can reliably estimate Id as: 

I D ~P D < D\cT^d\0 >< 0| dT^c\D > - 
4{P D ) a < D\ cYy, d |0 >< 0| dk a T^ c \D > +2 < D\ cT^k a d |0 >< 0| dk a T^ c\D > + 

2 <D|cr„d|0><0| dk 2 T^c\D> . (204) 

Let us define the matrix elements in Eq.(204) as: 

<0\qk v i'y 5 c\D(P)>=r D <x> P P u , (205) 

< 0\qK lfll5 c\D(P) >= i f D < x > A ^Pu-^g^uP 2 ^ , (206) 



<0\qk 2 llll5 c\D(P)>=tf D P^ <k 2 > A . (207) 

The quantity < x > has the meaning of the mean momentum fraction carried by the 
light quark (in the P z — > oo frame), and < k 2 > is the characteristic value of the light 
quark 4-momentum squared inside the D meson. So: 

^(l-3<.>,+f<^-2^|^). (208) 

Using the equations of motion for the matrix elements Eqs.(205), (206) it is not difficult 
to obtain: 



1/ (M c -m q ) 2 \ 4 \i^_n( A Q CD \ 

2 I 1 ' M 2 I <->^=3< r> " i + 0i — ] 



M c ) 



. (209) 
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As for the value of < k 2 >, the estimates obtained from the corresponding QCD sum 
rules for this quantity show that it is not far from its value for the vacuum quarks (see 
Eq.(42)). So, we have Q: 

<a;> A ^0.15, < -k 2 > A ~0AGeV 2 , <A^>~0.35M^. (210) 

It is seen that it is of great importance here to account for the spectator quark momenta, 
in spite of the fact that these are power corrections only in the formal limit M c —>■ oo. 
In what follows, we use the same estimate Eq.(210) also for the fig. 6b contribution: 



<D\ X 2 cT p — d-dT^— c \D>~ 0.35 M 2 D < D\cT p — d ■ dT p — c\D > .(211) 



— d-dT u —c\D >~ 0.35 MP, < D\cT u — d ■ dT „ — 
2 M 2 ' d \ n 2 >* 2 

The corresponding expressions for the B mesons look as: 

<x> A ~6%, < -k 2 >~QAGeV 2 , < A| >~ 0.8 M| , (212) 

so that the effect discussed is sizeable even for B mesons. 

13a . Calculation of the D ± '°' s decay widths 

Let us collect now the results obtained in previous sections (see also Appendix). 
As it was shown in sects. 3a, 8, keeping only the leading term and first corrections 
0(1 /M 2 ), we obtain a common nonleptonic decay width for all D ±,0,s mesons: 

C 2 M 5 
Y° nl ~ 1.59 r o , r o = -^— f- ~ 8.4 • 10- 13 GeV . (213) 

o4 TV 6 
Accounting for the contributions of the four-fermion operators, we obtained: 

I) The nonvalence contributions (which are the same for all D ±,0,s mesons (see 
sects. 3a, 10)) from the diagram of fig.4: 

5T PNV = ——— < D\ L PNV (fi ) \D >= 



g2 ^(i-v; 2/9 )m><d\ 



A-kM d 

G\ 



J a p (S)-c^ lp (N v + N al5 )c 



<l-V^ /9 )Ml 



itt } r L r- '•'" '"" 



\D > ~ 



--r 2 n (u )N v [ 1 i 

9 DK ^ J \ 2 Ml 



(214) 



*W _64 ^ 2 fl^Ml {i _ ^ ^ (,_M!)^ _ , 5 . (215) 



r o ~ 9 Ml v lco ' "\ 2M 2 D/ 
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That the effect is large for the D mesons can be seen from a simplest rough estimate: 



fci ~k 2 ~<x> A PD, X ^ (Pd - 2 < x > A Pd) ^ (1 - 2 < x >a) Md - 0-5 M 
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II) The nonvalence contribution from the diagram of fig.9 (see sect. 10): 
1 



5Tm V = 



<D\AL(p )\D> 



NV- 



(,/ " M 2 DPo n NV 



n 



NV 



2M D ' v ^ 0/l ' v AtxM d 27 

A u + A d (\V ud \ 2 + 6IK-I 2 ) - A a (t\V cs \ 2 + \V cd \ 2 



-9.3 



is = —~ —2 -0.64, 
where the subscript NV means nonvalence contributions, 



5T 



NV 



T 



~ 7T 

27 



32 , /B^o) M% 



Ml 



p VL 



NV 



-0.15. 



So, the total nonvalence contribution is: 



AIW _ [5T PNV ] + [5T 



NV 



•15%] + [-15%] = -30%, 



r r 

III) We have for the valence contributions (see sects. 3a, 9, 10 and figs. 6, 10): 



(216) 



(217) 



(218) 



AT nl {D°) 



1 



2M D 



< D°\AL^(p )\D° > v c -^{-O u M 2 D (C C G - P v ) , (219) 



4nM D 



^p.^l6n^O u (C' -n)c + 33% 



(220) 



^^ A- < fl +|AL< W+ >, 

J- o J- o 



16 7T 



2 f 2 D {Po)Ml 2 < A| > 



Mi 5 



2 M S 



\ 2 



48n'--^\V ud \ 2 O d [^^C c n -^P{ 



M„ 5 



M 2 ^ G 3 r 



-62%] + [ -21%] ~ -83% ,(221) 



where the subscript V means valence contributions. 
For the D s meson Q: 



P = p- < £> s |AL w (/x )|D s >, 



16 7T 2 |K 



J 2 MM% <X 2 D > 



Ml 



S,, 



M 2 D 



+ 



Mr 



16vr 2 |Vy 2 ^ O s (Q-P£ 



+ 



_48tt 2 |K s | 2 ^O 

I t»| M5 






V' 



-3%] + [-3%] + [-!%]- -7%. 



(222) 



16 1 am indebted to N.G. Uraltsev who pointed out an arithmetical mistake in the original calculation 
of this correction. 
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We have therefore for the nonleptonic widths (the experimental values are given in 
brackets, all values are given below in units of 1CT 13 GeV): 

T nl (D + ) -T [ 1.59 -0.30 -0.83] -0.46 r o - 3.9 {4.05}, (223) 

r nl (D°) ~T [1.59 -0.30 + 0.33] - 1.62 r o - 13.6 {13.5}, (224) 

T nl (D s ) - r o [ 1.59 - 0.30 - 0.07] - 1.22 T - 10.2 { ? } , (225) 
We have for the semileptonic widths (see sects. 4, 11): 

Y le P\D + ) = 1.06 (input) , (226) 

T izp\D°) - T lept (D + ) , (227) 

f lept (D s ) - r' ept ( J D+) [ 1 - 7% ] - 1.0 . (228) 
We have to add also to T lept (D s ) the "D s — ► r^" contribution: 



r lept '■ " M 5 



2 



AC; 

2 



M D 



! =» J6 /'«v 



,s J 



200 MeV 



T(D S -► r z/) - 0.4 , r^*(£> a ) - 2 • 1 + 0.4 - 2.4 . (229) 
So, we have finally for the total decay widths: 

r tot (£> + )~ [3.9 + 2.1] -6.0, {6.2} (230) 

T tot {D°)~ [13.6 + 2.1] -15.7, {15.6} (231) 

T tot {D s )~ [10.2 + 2.4] -12.6, {13.8}. (232) 

13b . Calculation of the B ±,0,s lifetime differences 

Substituting the corresponding numbers into the formula Eq.(25), we obtain: 

2CI + C 2 <B\bb\B> . . 

— + ~ - 1.127, ! ! -1.00, J rad -0.98, 233 

3 '2M S ' v ; 



Gl Ml\V r} v 

64 7T 3 



C-1.10r o z o , T = " F J\ cbl . (234) 
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The values of the phase space factors z* |39| (for M b = 5.04 GeV, M c = 1.65 GeV, M T 
1.784 GeV) are: 



C -0.460, zf~ 0.131, zf~ 0.105. (235) 

So, neglecting the 4-fermion operator contributions (which are small, see below) we 
have for the relative yields: 

n e J f ~ | (l) ud + (0.285) cs + 2 (0.260) ei/+M „ + (0.060)™ 1 - 1.87 , (236) 

and the branching fractions are respectively: 



Br (b — > ccs) ~ 15% , Br (b — > cev) ~ 13.!! 



70 



Br (b ^ctu)) -3.2%, Br (-) ~ 0.23 . (237) 

e 



The above value of Br (b — > cez/) is ~ 20% larger than the LEP data |26j , and this is 
a well known difficulty f|(| , [|4|] , [0 . What we can add is that the 4-fermion operator 
contributions calculated below are too small and can not help here. From our viewpoint, 
one of the important reasons for this discrepancy is that the radiative corrections to 
the nonleptonic widths are known for massless quarks only, and there are reasons to 
expect that accounting for nonzero final quark masses will increase considerably the 
nonleptonic radiative correction. Q Another evident reason which comes to mind is 
that the energy release is not large in the b — > ccs mode, so that it can be enhanced due 
to nonperturbative effects. Although the present data do not support this, it seems 
that the experimental numbers can change here with time. 

Let us proceed now to the calculation of the B meson lifetime differences. Using 
results obtained in preceeding sections (see sects. 3b, 9, 10) , we have for the valence 
contributions: 

AT(B") 2 <A|> 

■' 167T - 



S, 







3.2%, (238) 



M, 
M 2 B V'Ml 



£b = -" -f -^ ^ 1.31 . (239) 

< A| > A MS * ' 



1 



^5> 



17 To illustrate possible changes, we can put M c = in the leptonic radiative correction and obtain: 
Br (b — > cev) ~ 12.9% , in comparison with ~ 11.4% from the LEP data. 
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AF < B °>- ifa' (1 -,)>*0. 



r^j 



r v Ml 



:i + |)^l-^ 



0.6% , (240) 



AT(B S ) 1 „. „ U/2 M 



167r 2 (l-4a;) 1 / 2 ^|0 M [(l-x)C ( f-(l-2a;)P^J -0.5%, (241) 

where: x = M 2 /Mq . The nonvalence penguin contributions are also small, both 
factorizable (see fig. 4) and nonfactorizable (see fig. 9): 



AIW ^ 64 PM Ml ( Mj\ A ^ 

T 9 Ml \ 2Ml) 



5% . (242) 



Ar w 32 2 ft{n„)M% 



r^j 



r - 27^ mi poX 



< ss > 2 



4(1 - xfA d -A u - vT=4c (1 -2a;) — ,A U 

< uu > A 



-0.1%. (243) 



It is seen that, analogously to the D mesons, the largest effect is the negative contri- 
bution to the B~ width due to Pauli interference, but it is only — —3% here, while 
the B° and B s meson widths receive both only — 0.5% corrections. 

Let us comment finally on the accuracy of all the above calculations. It is extremely 
difficult to give a reliable estimate of the accuracy. Too many issues are involved, and 
estimates of various contributions vary in their accuracy. So, the accuracy of predic- 
tions for various quantities (the quark masses, \V c b\, the D meson decay widths, the 
B meson lifetime differences, etc.) differ from each other considerably. So, we do not 
even try to give here any concrete numbers, except for some "educated guesses" like: 

a) the c- and b-quark masses can hardly be less than 1.6 GeV and 5.0 GeV respectively; 

b) /b(M(,) can hardly be larger than 120 MeV; c) the value of \V c b\ can hardly deviate 
more than ±0.002 from 0.040 (with the experimental value of the semileptonic width 
taken at its central value, see Eq.(114)); d) the lifetime difference between the B~ and 
B° mesons can hardly be more than — 5%. 

14 . B° — B° mixing 

The effective Lagrangian which determines the mixing width, r 12 , is obtained 
directly from Eqs.(61),(65) by the substitution: (b...s) — > (s...b) Q: 



C 2 (V V*) 2 

L width\ M b) - T- J-^u L ]iu \ M b) 



L^\M b ) = S° u (sT,b)(rsT u b) + 0° u (s^-T,b)(s^-T u b), (244) 



18 It is implied in Eq.(244) that when calculating the matrix element each of two b-operators acts 
on both sides, and this gives the additional factor 2 which is compensated by the additional factor 
1/2 introduced into Eq.(244) 
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where q denotes the s or d quark field. Because the difference between A and A (see 
sect. 12) is not of great importance for the B mesons but complicates considerably 
the renormalization formulae, it will be neglected in what follows. Then, using the 
relation: 



A<\ 



A a 



s-T,b){7s-T u b)= ,, 



2 (sT fl b)(sT l/ b) + ^g^(sT p b)(sT p b) , (245) 



we can rewrite Eq.(244) in the form: 



L^\M b ) 



Q° n° 



sTfj,b- sT v b- -g pu sT p b- sT p b 



+ 



M h 



\{s°u + \o° u ) g»v[sT p b.sT p b\ 



M b ■ 



(246) 



The operators in the square brackets in Eq.(246) renormalize multiplicatively [[13 
that we obtain: 



so 



L^\^) = a [7sT p b-7sT v b}^+P[g pv sT p b-7sT p b] iio 



" a^ 3 (,v;;-|o° 



-1.93 



P 



4 
A bo 



A/. ( \o° u + si)+ aH 3 (Jo: - s° u 



0.83 



' <x.(M e y 



12/25 



a s (M c 



4/9 



1.614 . 



(247) 



Using now (see sects. 8, 9): 



<B°(p)\L^*\n )\B°{ P )> 



2 1 



1 



2/3-f 2 B (^)M B g pu , 
< B°(p)\ L^ x) ^o) \B°(p) > nonfactor ^ 2 



_2a^C b G 



P 



+ 



2g pv a (-C G + 2P^)+[3{-6C G + 8P } 



v 



(248) 



and contracting with: 



Tjff = Vl-4x|— ^ (p p p v - p 2 g pu ) +p 2 xg pi 






(249) 



we have: 



Iff < B°{p)\ L%*>Qio) \B°{p) > factor^ 2 (-0.57) f 2 M M B ~ -7.0 GeV e , 

T<& < B°{ P )\L^ x \n ) \B°(p) >nonfactor^ l.2GeV 6 . (250) 
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So, 


we 


have finally for 

£ KM 16tt 


ri 2 : 

T l2 {B) 




T 

2 f 2 B ^o)Ml 
Ml 



1 <■ R°\ r( mix ) I ro ^ 

1 <B ^ width Iff > ^ 

[ -0.57 (1 - 17%)] ~ (-2.7%) £ K m , (251) 

where £km is the ratio of the Kobayashi-Maskawa factors. It is seen that the nonfac- 
torizable contributions appear to be surprisingly large here and decrease the mixing by 
-17%. 

An analogous situation takes place for the mixing mass of the B° and B° mesons. 
The effective Lagrangian can be found in any review and has the form: 

L^(M b ) = C [sT p b-sT p b] Mb , (252) 

where C Q is a known coefficient (see sect. 15). The above operator renormalizes mul- 
tiplicatively p2"|, so that there appears only the factor A bo (see Eq.(247)) when it is 



renormalized to the point \x Q . The matrix element is given by Eq.(248), so that we 
have: 

< B°\ L^:}{M h ) \B° > factor ^ Abo \f B {lJLo)Ml C ~ 0.95 GeV A C D , 
< B°\ L^(M b ) \B° > n0 nfactor^ 2A 6o [-6C& + 8P^} C Q ~ -0.176 GeV 4 C Q .(253) 
On the whole: 

< B°\ LSS(M h ) \B° > M ~ A bo lf 2 B (Vo)M 2 B ( 1 - 18% ) C Q * 
^f B (M b )M 2 B (l-18%)C o -0.8GeV i C o ; B B (M b ) ~ (1 - 0.18) ~ 0.82 , (254) 
and the corrections to the factorization approximation are also very significant here. 
15 . The unitarity triangle 

The purpose of this section is to show that the results obtained above are marginally 
consistent with the available data and can be used to determine the parameters of the 
unitarity triangle (we use below in this section the Wolfenstein parametrization and 
the notations from |45|| ). 

We would like to emphasize that we have not tried in this section to account 
for all the various predictions available in the literature for the parameters involved 
(/b, B b , Bk, etc.). In fact, there are a number of papers which try to account with 
an equal weight for all the values available in the literature for these parameters. This 
results in large uncertainties which prevent to obtain more or less definite results from 
the available experimental data. Instead, we prefer to use only those results which, 
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from our viewpoint, are more reliable and this allows us to obtain the concrete predic- 
tions for the unitarity triangle parameters. 

1) . Using: V cb = A A 2 , A ~ 0.22 and (see Eq.(114)) V cb ~ 0.040 , one has: 

A -0.825. (255) 

2) . The Bd — Bd mass difference is given by: 



AM, 



G 



Xd = 



B 



6vr 2 



t b -\ M^ M B f B B B aj [nl B S{x* t )\ \V td \ 2 = (0.71 ± 0.07) , (256) 



M B 



S(x) = x 



9 



4 4(x-l) 2(x-l) i 



3 

+ 7: 



x 



2 Vx-1 



logx , 



x+ 



M 



M, 



') T„,r 2 = .4 2 A' 

w 



2 i „2 



[i-pY + v 



0.77-10" 



(l-p) 2 + V 2 -(257) 



Using in Eq.(256) ( see Eqs.(128), (254) and || ) 0: 



M t * ~ 180 GeU , [ff 2B S(x* t )] ~ 0.85 • 2.7 ~ 2.30 , 



(258) 



r B ~ 1.6 • 10 - " s , / B (M 6 ) ~ 113 MeU, B B (M b ) ~ 0.82 . 



(259) 



one obtains: 



2 i „2 



;i-pr+*j 



^rf 



rvj rs_; 



0.35 



2.0 



(260) 



3) . Using also: 



V ub 2 



V, 



cb 



A 2 (p 2 + r/ 2 ) ~ 1 • 10~ 2 , (p 2 + r/ 2 ) -0.2. 



(261) 



one obtains then from Eqs.(260), (261): 



p~-0A, 77 ~ 0.2, 8 = tg~ x (-) ~0.857r. 

P 



(262) 



sin (2 a) ~ 0.60, sin (2/3) ~ 0.28 , sin (2 7) ~ -0.80 . 



(263) 



19 In order to use the numericals from Ej| we use the mass M t * which is the so-called MS mass; 
M t * = 180 Gel^ corresponds to the pole mass: Mf° e ~ 190 GeT^ , which is marginally consistent with 
the CDF and LEP data, El. 
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4) . The CP- violating part of the K° — K° mixing can be written in the form [45| : 



e— / 4 e K ~ C e B K A 2 A 6 V \p o + A 2 A 4 (1 - p ) V * K S(x* t ] 



(264) 



a 



G 2 F F 2 K M K M 2 W 
6V2tt 2 AM k 



3.85 -10 4 



(265) 



Po = Xc ( V3 S 3 (x t ) ~ 77i ) , ^(Xt) 



. x t 3 X| fx t \ogx t N 

log — - 1 

x c 4 x t - 1 V x t - 1 



(266) 



■*■ /. 



M 



M, 



IV 



M c ~ 1.65 GW , M t * ~ 180 GeV 



(267) 



ry! ~ 0.80, 77*^ ~ 0.57, r] 3 ~ 0.36 . 
Substituting all this into Eq.(264), one has: 

e" i7r / 4 e K = 2.26 ■ 10" 3 ~ 7.3 • 10~ 3 B K 77 [ (1 - p) + u ' 



(268) 



(269) 



LOn 



\V cb \ 2 v*2 K S(x* t 



0.38 



77(1- p + w ) 



0.31 



B 



K 



(270) 



Let us recall g5j that, unlike B B (M b ) in Eq.(254), 5^ in Eq.(264) is defined as: 



<^°|sr^.sr,rf|K°> M =- J B^)/|.M| : 



(271) 



B K = B K (p)a; 2 ^(p). 



(272) 



The characteristic value of Bx obtained from lattice calculations is: Bk = (0.9 ± 0.1) 
p7| , and substituting B K ~ 0.9 into Eq.(270), one obtains: 



7] (1 - p + cj ) ~ 0.345 



(273) 



Using now: p ~ —0.40, uo ~ 0.38 one obtains 77 ~ 0.19, in agreement with Eq.(262). 
5) . Clearly, we can proceed in the opposite way: using B K ~ 0.9 we obtain from 
ex and Xd- 



77 (1 - p + 0.38 ) ~ 0.345 , (1 - pY + 7] 2 ~ 2.0 . 



(274) 
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This gives us then: p ~ —0.4 , 77 ~ 0.2 and | V^/Vy ~ 0.10 . 

With the above parameters, the CP- violating asymmetry in B° — > \P K$ decay is: 



A(B°->VK S ) 



*d 



1 + x 2 , 



sin(2/3) -0.13. 



6) . The quantity e'/e can be represented in the form [f48f : 
Im\ t 



€ 

e 



1.7 



p(l/2) _ p(3/2) 



0.4-10" 



p(l/2) _ p(3/2) 



(275) 



2 \5 



7mA 4 = 77/r A ~0.7- 10 



(276) 



where P^ 1 / 2 ) and p( 3 / 2 ) are expressed through the known coefficients di(M t ) and the 



parameters P, 



(1/2) „(3/2) 



(. 



Bk 



These latter are determined from the corresponding ma- 



trix elements of the penguin operators, see [48|. In the factorization approximation: 



p(V2) 



B, 



(3/2) 



1, and the values obtained from both the lattice calculations and the 



1/N C expansion agree with the factorization approximation within ±20%. Using the 



values of di(LO) from the Tables 14, 15 in ESI and B { 



,(1/2) 



B, 



(3/2) 



1, one obtains: 



p(l/2) ^ 5>2 p(3/2) 



6.2, -~ 
e 



-0.4-10 



-i 



(277) 



The above value can not be taken too literally because it is a result of strong can- 
celations. To be conservative, the realistic value of e'/e can be considered to lie 
in the interval: ~ ±(1 — 2) • 10~ 4 . In any case, this is inconsistent with the NA31- 
result: e'/e = (23 ± 7) ■ 10~ 4 , while there is no contradiction with the E731-result: 
e'/e=(7±6)-10- 4 , 0,0,11. 



7) . The box diagram contribution to (Kl — K$) the mass difference, AMk 



M(K L ) - M(K S ), is |H, [H: 



{AM K } box = pLf K M K Ml A 2 B K Vl 

07T Z 



O(10 



-2\ 



(278) 



where the correction O(10 2 ) represents the t-quark contribution and other small cor- 
rections. Using: f K ~ 162 MeV, M c ~ 1.65 GeV, B K ~ 0.9, one obtains: 



{AM K } box ~ (4.0 • 10' l5 GeV) B kVi 
3.6 • 10~ 15 GeV) m ~ 2.9 • 10~ 15 GeV , (279) 

for 771 ~ 0.8. The experimental value is: AMk = 3.5 • 10~ 15 GeV, so that there is not 
much room for additional large distance contributions. 



20 These are due to the large value of the t-quark mass which enhances the electroweak penguin 
contribution p( 3 / 2) , see El . 
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16 . Summary and conclusions 

Let us summarize first in short the results obtained for the D mesons. 

It is seen that the overall picture is sufficiently complicated: there are a number of 
contributions, all of the same order and of different signs. As to the four-fermion op- 
erator contributions we were mainly interested in, the qualitative picture is as follows. 

Although they are formally 0(1 /M%) corrections, these contributions are very im- 
portant numerically and comparable with the Born term. In particular, their final 
effect in the nonleptonic widths is much larger than those of the leading 0(1 /M 2 ) 
corrections. There are two reasons for this. 

i) Various 0(1/M 2 ) corrections, being only a few times smaller than the Born term, 
cancel each other strongly, see Eq.(147); 

ii) The four-fermion operator contributions are the first to gain a large numerical 
factor from the larger two-particle phase space. Otherwise they would have been much 
smaller than separate 0(1/M 2 ) terms. It is clear that this effect operates one time only. 
So, there are reasons to expect that all other 0(1/M 3 ) and higher order corrections 
are much smaller, and those considered in this paper are the main ones. Some simple 
estimates confirm this. 

As for the relative significance of various contributions, the picture is as follows. 

1) Most significant is the destructive Pauli interference effect (see fig. 3c) of two 
d-quarks in D~ decay (~ —60% of the Born term, T ). 

2) The (cross) annihilation contribution (see fig.6b) ensured by the nonperturbative 
nonfactorizable gluon interaction decreases further (~ —20% of r o ) the D~ nonleptonic 
width. 

3) Both the above contributions could have appeared much larger, but are strongly 
suppressed by the much smaller two particle phase space. The reason is that the 
"normal" total 4-momentum of the quark pair is (see fig. 3a): A ~ (P c + k\) = Pd (k is 
the momentum of the spectator quark). However, it is: A ~ (P c — A^) = (Pd — k± — k?) 
for both these " cross-contributions" . Because the charm quark is not very heavy and 
spectator quarks carry a significant fraction of the charmed meson momentum, this 
leads to a strong suppression: A 2 ^> A 2 . This effect remains noticeable even for the B 
mesons. 

4) The (direct) weak annihilation contribution (see fig. 6a) increases significantly 
(~ 30% of r o ) the D° nonleptonic width 0. 

5) The nonvalence penguin contributions, both factorizable (see fig. 4) and non- 
factorizable (see fig.9)), are very significant and, on the whole, diminish considerably 
( ~ —(20 — 25)% ) both semileptonic and nonleptonic widths. 

6) There are no noticeable positive valence contributions into the D s meson width. 
Those which are available are negative and decrease (~ —7% of the Born terms) both 
semileptonic and nonleptonic widths. As a result, there is a sizeable difference (~ 10%) 
between the above calculated and the experimental numbers. A possible explanation 



21 Let us emphasize that the annihilation contribution ~ 30% into the inclusive width does not 
contradict that it can be ~ 100% in separate exclusive modes, see p5[ . 
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may be due to the SU(3) symmetry breaking effects which were neglected in these 
calculations. At first sight, however, most of them tend to increase the discrepancy 
rather than to decrease it. Clearly, this subject requires careful investigation which is 
out the scope of this paper. 

One of the important results of all the calculations above is that the factorization 
approximation works well for the matrix elements, i.e. the nonfactorizable contribu- 
tions which have the same parametrical behaviour at large Mq are, in comparison, an 
order of magnitude smaller. Nevertheles, they are of importance due to specific fea- 
tures of the problem considered. Firstly, these nonfactorizable parts enter the effective 
Lagrangian with much larger coefficients. Secondly, the factorizable contributions are 
additionally suppressed by the smaller phase space (see point 3 above). 

Let us add a few words about baryon lifetimes. It seems clear that the pattern 
here is very unlike those for the pseudoscalar mesons. The entire structure of the 
matrix elements of L e ff is quite different, and the reasons connected with the helicity 
suppression of factorizable contributions are not operative here. So, one can expect 
that: 1) the scale of the 4-fermion operator contributions to the baryon lifetimes is 
potentially a few times larger (term by term) than those for the pseudoscalar mesons; 2) 
the significance of nonfactorizable contributions will be smaller for baryons. Therefore, 
it seems, the main problem here will be to calculate reliably the factorized parts of the 
matrix elements. 

For the B mesons, all the above qualitative properties remain true but , clearly, 
the role of all power corrections to the Born term becomes much smaller. It seems 
clear now also that the factorization approximation will be especially good for the B c 
mesons. 

It is important that we understand now the properties of all the main contributions 
giving rise the lifetime differences and that there is sufficiently good agreement between 
the calculations for the D mesons and the experimental data. This gives us confidence 
that the predictions obtained above for the B mesons are reliable. Therefore, we can 
insist now that the lifetime difference between the B^ and B° mesons will not exceed 
~ 5% , while those between the B° and B s mesons have to be even smaller. As it was 
pointed out above, the scale of the 4-fermion operator contributions for the b-baryons 
can be a few times larger. 

Surprisingly large corrections to the factorization approximation (~ —18%) are 
found for B° — B° mixing. This reduction of the B^ — Bd mixing mass is of importance, 
in particular, for the parameters of the unitarity triangle. These latter are calculated 
in sect. 15 and are: 

A ~ 0.22, A -0.825, p ~ -0.40 , r? ~ 0.20 , (280) 

sin(2a) ~ 0.60 , sin(2/3) ~ 0.28 , sin(2 7 ) ~ -0.80 . (281) 

We would like to emphasize finally that the experimental value of Br(D — > ev + X) 
requires a noticeably larger value of the c-quark mass (M c ~ 1.65 GeV), in comparison 
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with those (M c ~ 1.4 — 1.5 GeV) in common use. As a result, because the mass 
formulae tell us that the quark mass difference is close to those of the mesons, this 
leads to a value of the b-quark mass (M?, ~ 5.04 GeV) which is also considerably larger 
than those in common use (Mb ~ 4.6 — 4.8 GeV). Further, the chain of the argument 
is as follows. The large value of Mb leads to a small value for fs '■ J'b(Mb) — 113 MeV. 
Together with (see sect. 14) Bb — 0.82, this requires a large value of the t-quark MS- 
mass: M 4 * ~ 180 GeV, to reproduce the experimental value of Xd — 0.70. This results 
finally in the rather small value of e'/e ~ (±1 ■ 10~ 4 ), which is incompatable with the 
NA31-group result but does not contradict to those from the E731-group. 
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Appendix 

For the convenience of the reader, we list below the notations and numerical values 
of various parameters used in the text. 

Parameters: 



M b = 5.04 GeV [Eq. (103)] , M c = 1.65 GeV [sect A] , /i 2 , = 0.5 GeV 2 [Eq. (43^82) 
a s (M%) =0.204, a s (M 2 ) =0.310, a s (fi 2 ) = 0.580 [Eq.(49)]', (283) 



M 2 
<0\qi l5 Q\D>=r D = f D1 ±; (284) 



f D (fi ) = lUMeV, f D (M c ) = l65MeV [Eq.(123)\; (285) 

f B (fi ) = 89 MeV, f B (M b ) = 113 MeV [Eq.{12S)); (286) 

<Xs(l£) 



ri 



CO 



a s (Ml 



1.87, Tco = rj£ = 1.37; (287) 



^^-284 n*-^? 
a s (M 6 2 )-'- M ' Vbc - a s {M 2 B) 



Vbo = "'^ = 2.84 , ^ c = " a .™ =1.52, r 6o = % 2 C 7/50 r co = 1.715; (288) 



1 -ik 2 / 9 ) =0.13, ^ = ^1 Jog 4= = 0.1; (289) 



id in 



^-^ <J " 8t <^%'? e|l?> -S jffl ^-^ <-» 

JV„ = 2.54 [£gs.(51),(52)], A = 0.40 [£gs.(73), (74)] . (292) 

The matrix elements (sects. 9, 10): 
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2 
< D(p)\cT fM q-qT u c\D(p) > Mo = p^p v ff,([i ) 5 V + ^p g^ur 2 D (fi ) , (293) 



< D{p)\cY tl —q-q—Y v c\D{p) > Mo = f -^ ~ 9^\ C G 5 V + g^Py 5 V - 

■^Po9puTd(i*o), (294) 

Here: Sy is unity if the quark field in the matrix element has a valence flavour and is 
zero otherwise. 

C C G = 0.48 • f (T 2 GeV 4 [Eq.{162)}, C b G = 1.15 • 1(T 2 GeV 4 [Eq.(163)] ; (295) 
P^ = 0.685 -lO^GeV 4 , P v = 0.18 • 10~ 2 GeV 4 [Eq.(192)\; (296) 



A a ,_,, 4 , , , A M 2 

r2 



< ^(/')! (•/;; • r-,„(] dz 7 5) T c| J D(p) > Mo = --r 2 D (p ) I 1 - ^ ] + • • • . (297) 



where the dots mean that we neglected corrections because this matrix element appears 
itself as a correction in the above calculations. 

In addition to the 4-quark operators considered in the text, there is also the nonva- 
lence 4-quark operator contribution originating from the Born diagram, and we show 
below that it is negligibly small. Using the Fock-Schwinger gauge: 



1 1 

-X v G vli (0) + -: 



x^A^x) = 0, A^x) = -x v G vli (0) + -x u x x D x G^(0) + ■■■ , (298) 



the Born diagram contribution can be represented in the form: 

T B ^ n ~<D{p)\{c{x)*d*{l- lb )-dl{id)c{x)} x jD{p) > . (299) 

Using the equations of motion: idc = M c c — g s A c, it is not difficult to obtain: 
T B om ~ Ml < £>(p)|c(0)(l - 7 s) (1 - ig s ^uG^{0)) c(0)\D(p) > - 
l -gl Ml < D(p)|c(0) 7m (1 + 75) yc(0) • J a ,(0)\D(p) > + l -8 v + l -8 A + ■ ■ ■ , (300) 

<— » — > i — 

where " • • •" denotes higher dimension terms. Also (8^=8^ — d^), 

8 V =< D{p)\c(0)^duc(0) ■ GlM\D{p) >= (301) 



51 



due to C-parity, and 



X a 



5 A =< £>(p)|c(0)7 M75y d*c(0) ■ GI„(0)\D(P) >= (302) 

due to P-parity. Therefore, the correction due to the 4-fermion operator is: 

*4 = 1 - ^ra s (M c ) 2M ^ M3 < £>|c(0) 7/ .(l + 7 5 )yc(0) • J£(0)|D > . (303) 
This becomes after the factorization of the matrix element: 



4 fu^ jjW /, M c 



rt'j ~ I + -^n,(..U r" ;;,' | 1 - ^J * 1 + 1-3 • 10"*, (304) 



so that it is negligble even for the D-meson. 

For the semileptonic decays the above described 4-quark operator contribution is 
the only one, in addition to the cc, co^yG^yC and 4-quark operators described in the 
text, which is capable to give ~ 0(Aq CD /M^) correction. For the nonleptonic decays, 
in addition, it is not difficult to see that it is sufficient to replace A^ in the last term 
in Eq.(25) (which is due to the fig. 2 diagram) by its original expression: 

<D\c% W id v c-G%\D> 
Ag "> 2 < D\cc\D > W • (3 ° 5) 



52 



Figure captions 

Fig.l The Born contribution 

Fig. 2 The soft gluon emission giving rise the correction 0(1/Mq) 

Fig. 3 The diagrams contributing into the 4-fermion operators 

Fig. 4 The nonvalence factorizable penguin contribution (PNV) 

Fig. 5 The diagrams contributing to the sum rule Eq.(118) 

Fig. 6a The nonfactorizable gluon contribution to the weak annihilation 

Fig. 6b The same for the cross weak annihilation 

Fig. 7 The diagrams for the meson transition into a current 
with an emission of a nonperturbative gluon 

Fig. 8 The nonfactorizable contribution to the matrix element < D\S 11V \D > 

Fig. 9 The nonvalence nonfactorizable penguin contribution (NV) 

Fig. 10 The valence nonfactorizable penguin contribution 
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